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Verallgemeinerung und Vereinheitlichung der Wirbelsätze 


ebener und rotationssymmetrischer Flüssigkeitsbewegungen 
Von ©. Truesdell, Universität Maryland, USA. 


Es wird ein einfacher, rein kinematischer Erhaltungssatz des Wirbelbetrages bewiesen, der für Bewe- 
gungen gilt, die „komplex-lamellar‘‘ sind und noch eine weitere Bedingung erfüllen. Dieser Satz umfaßt 
die’aus der klassischen Hydrodynamik und der Gasdynamik bekannten Wirbelsätze für ebene und rotations- 
symmetrische Flüssigkeitsbewegungen und liefert eine Verallgemeinerung für sie. 


A simple kinematical conservation theorem for a scalar magnitude associated with the vortictiy is shown 
to be valid for complex-lamellar motions satisfying a certain further condition. This theorem includes and 
generalizes the known vortex theorems for plane and rotationally symmetric motions, both for classical hydro- 

ee and for gas dynamics. 


On d&montre un simple th&or&me ein&matique, d’apres lequel une quantite scalaire liee avec le vecteur 
tourbillon se conserve dans tout mouvement complew-lamellaire qui satisfait aussi a une certaine condition. 
Ce iheor&me renferme et generalise les theor&mes tourbillonaires dejü connus pour les mouvements plans 
ou jouani de symetrie de revolution, et dans U’hydrodynamique classique et dans la dynamique des gaz. 


okassIBaeTca IPOCTaA UNCTO KEHeMATHYECKAA TEOPeMa O COXPAHCHHH BEIIMYNUHEI BUXPA 
CKOPOCTH. MeÜCTBUTeNBHAA MIA ABEKeHHH, IIUXMHH BUXpeli KOTOPBIX IEePHEeHAHKRYIAPHEI K 
JHHHAM TOKA MH YMOBAIETBOPAMIMUX Eile OMHOMy AaulbHefilemy yCAIOBNW. JTa TEeOpeMma 


3 BEIIWYAET B CEÖH HU O60ÖIMAET HSBECTAHE TEOPEMEI KIACCHYeCKOÄ TUNPONHHAMHKH U TA30BOH 
# AHHAMUKRH AI INIOCKUX H OCECHMMETPHYHEIX ABIMKeHÄ 3EHNKOCTL. 

= 5 : 1. Die bekannten Wirbelsätze 
F -  Sehon J. L. d’Alembert, dessen Arbeiten große Fortschritte für die analytische Hydro- 
—  dynamik einleiteten, hat zum ersten Maleinen Wirbelsatz entdeckt. In einer Auseinandersetzung 


mit L. Euler über die Möglichkeit von Wirbelbewegungen*) bewies er folgende, hier etwas 
allgemeiner und moderner formulierte Behauptung ?)?): 


D’Alembertscher Wirbelsatz [1761]: 
Gegeben sei eine stetige, ebene Flüssigkeitsbewegung mit der Geschwindigkeit d, der Be- 
 schleunigung a und der Dichte 0. Es sei 


ED To > a erh 


der Wirbelvektor, wsein Betrag. Der Beschleunigungsvektor a erfülle die d’Alembert-Euler- 
sche Bedingung ®) 
En EN 1 ER EN EEE 7 y 


Dann gilt während der Bewegung für jedes Flüssigkeitsteilchen: 
Ba 2 RN RI DEE RE IP Pe a RE “(3% 
0 


Ganz analog lautet der 


1) Die Bezeichnung ‚„‚Wirbelbewegungen‘‘ wurde erst viel später von Stokes (1845) und Helm - 

holtz (1858) eingeführt. 
? 2) Die Literaturhinweise werden durch die Jahreszahlen der betreffenden Veröffentlichungen gegeben. 

- Man vergleiche hierzu das Verzeichnis am Ende der Abhandlung. ; u 
Be: 3) [1761, $XIII]. Der d’Alembertsche Beweis gilt nur für die Bewegung einer inkompressiblen 
_ Flüssigkeit, deren Geschwindigkeitsvektor die Form v(t,£) = f(t): v(t) hat. In den Lehrbüchern wird dieser 
Satz Helmholtz zugeschrieben. | 
Be Rus‘ 4) Sowohl J. L.d’Alembert [1752, $ 86], [1761, $8$ X, XV] (im Falle der ebenen Bewegung) als auch 
-L. Euler [1757, $XXXV], [1770, 88 42—43], [1806, $$ 154—155] (im allgemeinen räumlichen Fall) gründen 
4 ihre Theorie der Flüssigkeitsbewegung auf dieser fundamentalen Annahme (2). Nach E. Beltrami [187l, 
812] ist diese Bedingung notwendig und hinreichend dafür, daß sich die (schon früher von Hankel[1861] 

"und Lord Kelvin [1869] eingeführte) Zirkulation jeder geschlossenen flüssigen Linie während der Bewegung 
it ändert. Unter der Voraussetzung konservativer.Massenkräfte ist die Bedingung erfüllt, wenn die Flüssig- 
t ohne Reibung betrachtet werden kann und eine Gleichung der Form f (p, oe) = 0 gilt oder wenn sie inkom- 
ssibel ist, einen konstanten Zähigkeitskoeffizienten besitzt und überall die Gleichung rot rot w = 0 erfüllt 


nscher bzw. Craigscher Satz). 
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Svanbergsche Wirbelsatz [1841P): 
Gegeben sei eine stetige, rotalionssymmetrische Flüssigkeitsbewegung um die z-Achse, 
die die d’Alembert-Eulersche Bedingung (2) erfüllen möge. r sei der Abstand von der z-Achse. 
Dann gilt während der Bewegung für jedes Flüssigkeitsteilchen 


w 
= CONSEN 022,9 st 2 Ve Re) Ru LE 
Fi cons (4) 


In den eleganten Helmholtzschen Wirbelsätzen finden bekanntlich die beiden obigen 
Ergebnisse ihre räumliche Verallgemeinerung. 

Folgende Frage ist jedoch bisher noch nicht beantwortet worden: N 

Frage 1: Gibt es eine Klasse von Flüssigkeitsbewegungen, die die ebenen und rotations- 
symmetrischen, aber nicht notwendig alle räumlichen Fälle umfaßt und für die ein finiter, 
skalarer Wirbel-Erhaltungssatz von der Form “ 


mE een) 
(h = geeignete skalare Funktion) gilt? — — { 

Bei wirbelhaften Gasbewegungen ist die d’Alembert-Eulersche Bedingung (2) im allge- 
meinen nicht erfüllt; deshalb gelten die Helmholtzschen Sätze im allgemeinen nicht. Jedoch 
hat L.Crocco®) in einem Spezialfall zwei interessante, den rein kinematischen Sätzen von 
d’Alembert und Svanberg formal ganz analoge dynamische Wirbelsätze entdeckt. Diese 
wurden vor drei Jahren von R. Prim?) wesentlich verallgemeinert und lauten: 


1. Grocco-Primscher Satz [1936, 1947] 


Gegeben sei ein ideales Gas (p=o- RT) mit konstanter spezifischer Wärme in einer 
stetigen, ebenen, stationären, massenkräftefreien Bewegung ohne Wärmeleitung. Es sei v, die 
unter diesen Voraussetzungen vorhandene und längs jeder Stromlinie konstant bleibende Grenz- 
geschwindigkeit; es sei d, der sogenannte Croccosche Vektor 


Ve N Bar re 


sein Wirbelvektor sei 
b 10, > 108 Du se Bl N een RENT RR 


p sei der Druck. Dann gilt längs jeder Stromlinie: 
We n 
a SONST HE DEN TERN RB ER ES)" 
: @) 


2. Crocco-Primscher Satz [1936, 1947] 
Es sei eine um die z-Achse rotationssymmetrische Gasbewegung gegeben, die die entspre- 


chenden Voraussetzungen des ersten Satzes erfüllen möge. Dann gilt längs jeder Stromlinie: € 
W, je. 
= | RR N FE ERS ER a. Net 
AR cons (9) 
Somit ergibt sich eine der Frage 1 ganz analoge, bisher noch nicht beantwortete | 4 


Frage 2: Gibt es einen finiten, skalaren Wirbelerhaltungssatz von der Form 


w 

ET Eh A E10) 

für eine die ebenen und rotationssymmetrischen Fälle einschließende Bewegungsklasse? 
Außerdem ergeben sich zwei weitere Fragen. A 


KIT 


En Hi erst später von Emm ons [1944, App. I] und Vazsonyi [1945, $$ 
) [1947]. | 5. EI Der Dar A 
SE F REST 


RE) De ee a 
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Frage 3: Gibt es, einen den zwei Crocco-Primschen Sätzen bzw. der Gleichung (10) 
analogen Erhaltungssatz von der Form 


We 
hF(p) 


der bei Bewegungen nicht idealer Gase gültig ist ? 


OBEN ee NETTE Re (LL); 


‚Frage 4: Gibt es einen einfachen, direkten Zusammenhang zwischen den kinematischen 
Wirbelsätzen von d’Alembert und Svanberg und den beiden dynamischen Wirbelsätzen von 
Crocco und Prim? 

In diesem Aufsatz beantworte ich diese vier Fragen im positiven Sinne durch einen einzigen, 
neuen, rein kinematischen Satz über den Betrag des Wirbelvektors i,. Dieser schließt (5) 
und (11)als Spezialfälle ein und liefert somit eine doppelte Erweiterung der bekannten Wirbel- 
sätze: 

1. Eine Erweiterung der Bewegungsklasse. 

2. Eine Erweiterung der physikalischen Flüssigkeitsart, die eines skalaren Wirbel- Erhal- 
tungssatzes fähig ist. 


II. Der neue kinematische Wirbelsatz 


Seit H. v. Helmholtz und Lord Kelvin bildet man häufig Wirbelsätze mit Hilfe von 
Flächen- und Linienintegralen. Da jedoch ein finiter Wirbel- Erhaltungssatz für Flüssigkeits- 


teilchen, also mathematisch für einzelne Punkte gesucht wird, gehe ich hier von der Aufstellung 


einer Differentialgleichung aus. Dies entspricht dem Verfahren, das J. L. Lagrange®) bei der 
Entwicklung des Theorems vom Geschwindigkeits-Potential benutzte. 

Um den neuen Wirbelsatz in einer Form herzuleiten, der ihn für mehrere ganz verschiedene 
dynamische Anwendungen geeignet macht, seien zunächst einige mathematisch bequeme Hilfs- 
Parameter eingeführt, die keine physikalische Bedeutung haben. 

Es sei v, eine für jedes Flüssigkeitsteilchen konstant bleibende Funktion: 


Di, 


? 0 (12) 
Es sei 
v - 
d, En ? . (13) 
der allgemeine ‚‚Croccosche Vektor“, 
TR OL a ER Ba ne eh are BU far SR LAN 


der ,„Groccosche Wirbelvektor“. A sei eine skalare Funktion. Hain läßt sich durch einfache 
vektorielle Transformationen zeigen, daß folgende Gleichung gilt: 


1 or, ee, 
rot 4 ‚tm xD — grad = ar 5 TH ITER WE Re. 
+0.:.7 (Am) — Aw. 0.4 Atdivv,— Er 
Sei nun A als Lösung der Gleichung 


a) 0, 1m, 9A 
_ TE en lane e ALO 
m FE 3 ot vo. 9 126) 
gewählt, und sei » +0 irgend eine Lösung der Gleichung?) 
43 Dlog Pi 5 
; RD: div v -. (17) 
Dann folgt aus (15)10) 11) 
| Kar (Eon, 2)! = % a )- RN ER UR e 
„rot |d er + m. X »,Dt 5 5 vv = div(Amw.) - - (18). 


8) (1781, 8$17—19.] 
®) Von solchen Lösungen gibt es für jede stetige Bewegung unendlich viele. 
_ 10) Es sei bemerkt, daß der Spezialfall A = v, = 1 die grundlegende kinematische Wirbeldifferential- 
jung von E. Beltrami [1871, $6] liefert. Die sogenannte „Helmholtzsche‘ Gleichung, die daraus 
Bikpetzung von ‚Gleichung (2) entsteht, wurde schon von J.L. Lagrange [1760, Ch. XLII] ent- 


1 Ich danke Herrn Diplom-Fhysiker E.A. Müller, ‚ Göttingen, für seine Kritik an dieser Stelle. 
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Um von hier aus zu einer skalaren Gleichung zu kommen, bilde ich das Skalarprodukt der Gl. (18) 


. Aw 
mit —— und erhalte12); 
y 


A D 
4 = 1 "pp .): en -H+ Ele (vd. 1.) (mw, grad A) 


DIN 
Am: 410% a 
+: rot em x oo) | 


0 


(19). 


Diese neue kinematische Gleichung ist der Schlüssel zur Beantwortung aller 
oben gestellten Fragen und zur Vereinheitlichung der Aussagen über den Betrag 
des Wirbelvektors. 

Von hier an setze ich voraus, daß das Geschwindigkeitsfeld ‚‚komplex-lamellar‘* ist. Dies 
bedeutet, daß es Normalflächen besitzt. Nach Lord Kelvin!3) muß für solche Felder die 
Gleichung 

mV = 0, en nl re 


gelten. Der Wirbelvektor steht also senkrecht auf der Stromlinie. Bekanntlich sind die ebenen 
und rotationssymmetrischen Vektorfelder Spezialfälle dieser Klasse. Es folgt auch unmittelbar, 
daß die Normalflächen des Geschwindigkeitsfeldes Wirbelflächen und umgekehrt die Strom- 
flächen Normalflächen des Wirbelvektorfeldes sind. Das Wirbelvektorfeld mw ist also ebenfalls 
komplex-lamellar. Mit 
w=rot (2,8) = 910, + grad X obs ae ein 
folgt die reine Identität 
1 BR RE Er 


Also ist auch der Groccosche Vektor einer komplex-lamellaren Bewegung komplex-lamellar. 
Man sieht außerdem, daß auch die Umkehrung dieses Satzes richtig ist. 
Ich setze jetzt die Gültigkeit der Gleichung 


| [2 1 od, ) Ei | 
iv, rot (4 Free + 1m. X v & I A ea EEE 
voraus. Dann geht Gleichung (19) über in 
1%D (det F iv, Aw, | 
——|—)!=[——.YV ER RE Fr RE Er x AHLDERN 
20, Di\ » SEHR v) v a 
Umgeformt ist dies 
D Aw, m & ) 
a art EI KB IT OL 25 
% Dı°8 v We \Wg Er (25) 


Dabei ist die rechte Seite von (25) diejenige Komponente des in Richtung von iv, genommenen 


0 
Richtungsdifferentialquotienten ze, die in Richtung von iv, fällt. Ich wähle nun als örtliche 


krummlinige Koordinaten die Wirbellinien (die Koordinaten auf ihnen seien mit x, bezeichnet) 
sowie zwei zuihnen und untereinander orthogonale Kurvenscharen mit den Koordinaten % und %;. 
Dann ergibt sich für das Linienelement 


de—=hdatMmdathde.. 2.222222. (6) 
Benutzt man diese Definitionen in Gleichung (25), so erhält man nach wohlbekannten Formeln 


AN 
1D A RAR, _% dh .m0R 


era = ———. 427% 
VD on Da ÄIhro, 00,  Ihev, 0m 122) 
Mit den Gleichungen für die Geschwindigkeitskomponenten & 
Dx Dx DE 
»—=h7,; v=h7,; m Bi Eu a (28) 


folgt somit aus (27): 


DI: Au 00 \v lee "oh ologh | 
De Eh, Nm, Vi 


12) Die Anregung hierzu stammt aus einer speziellen und nicht ganz befriedigenden Arbeit von J. War - 
ren [1870]. Die Untersuchungen von Warren sind für nicht ebene Bewegungen fehlerhaft. 
18) [1851, $$ 69, 75.] 


Cut day, also 
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Macht man noch die einschränkende Annahne 
| oh 
ot 


(sie bedeutet, daß die Wirbellinien stationär sind oder, was für ein komplex-lamellares Geschwin- 
digkeitsfeld dasselbe ist, daß die den Wirbellinien normalen Stromflächen stationär sind), dann 
verschwindet das dritte Glied rechts in (29), und wegen (20) verschwinden auch die beiden 
anderen. Man erhält also einfach 


Re... Re er. (303 


Aw 
ons nen. en... ah 


für jedes Flüssigkeitsteilchen. 
Die Zusammenfassung der, vorangehenden Untersuchungen ist 


Der neue kinematische Satz über den Betrag des Croccoschen Wirbelvektors. 


Gegeben sei eine stetige Flüssigkeitsbewegung, die in stationären, den Wirbellinien nor- 
malen Stromflächen x, = const. stattfinde (d.h. eine komplex-lamellare Bewegung mit statio- 
nären Wirbellinien); es sei A der normale Abstand zweier benachbarter Stromflächen x, und 


d2—= Mday + '"* a A (W 1), 
es sei v, eine beliebige, substantiell konstante Funktion: 
Do, 
Di 0 (W 2) 
Es seien d, und iv, durch die Gleichungen 
v 
u =e—, VEETDU DE Ne enteo 6 teure (WAR) 
%o 
definiert; es sei A eine Lösung der Gleichung 
= O0, 1mW.09A 
d(— ae I ee 
sta (@ , ot % 0 Sr 
und » eine Lösung der Gleichung 
1 0 j 
TE RA real BE SM 8 PAS ea NEE Le: GE ANDEN 


Wenn es dann möglich ist, v, und A so zu wählen, daß sie außer (W 2) bzw. (W4) noch die 
Gleichung 
we rot [a2 Det m x 2.) = 0 BEN LER) 
U, 0 Re E 
erfüllen, dann ist während der Bewegung die Funktion 
Aw, 
hv 


v 


für jedes Flüssigkeitsteilchen konstant. 


Es sei bemerkt, daß nur die Bedingung (W 6) einschränkend ist, da man stets unendlich 
viele Funktionen, die die übrigen, nur definierenden Bedingungen befriedigen, finden kann. 


III. Erste Anwendung: Ein neuer Wirbelsatz der klassischen Hydrodynamik 


Es sei die d’Alembert-Eulersche Bedingung (2), die die klassische Hydrodynamik 
charakterisiert, erfüllt. Nach J.L.Lagrange!#) gilt die rein kinematische Identität 


CE v2 
BE tuxrdteadn 22222222222.) 
daher ist die Gleichung 
’ ef) - 
, rot (+ wxn)=0 a ER ER RFERFERSTEH |: 


der d’ cher Eulerschen Bedingung (2) äquivalent. Die Bedingungen (W2) und (W4) 
werden durch die Wahl w=1, A=1 befriedigt; damit folgt %,=d, ww, = 1; somit ist 


“) (1781, $14.] , 


bis; a 
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wegen (33) auch (W6) erfüllt. Setzt man » = 9, so wird auch (W5) befriedigt. Dann folgt auf 
Grund unseres kinematischen Wirbelsatzes 


Der neue Wirbelsatz der klassischen Hydrodynamik. 


Gegeben sei eine stetige Flüssigkeitsbewegung, die in stationären, auf den Wirbellinien senkrecht 
stehenden Stromflächen stattfinden und die d’Alembert-Eulersche Bedingung (2) erfüllen 
möge. Dann ist während der Bewegung für jedes Flüssigkeitsteilchen 


wobei die Funktiön h durch (W 1) definiert ist. j 

Im Spezialfall einer ebenen Bewegung folgt nach (W1) A=1; daher ergibt (34) den 
d’Alembertschen Wirbelsatz. Im Spezialfall einer rotationssymmetrischen Bewegung folgt 
gemäß (W 1) A=r, Gleichung (34) ergibt also den Svanbergschen Wirbelsatz. 

Damit ist Frage 1 beantwortet. 


IV. Zweite Anwendung: Ein neuer Wirbelsatz der Gasdynamik 
Ich betrachte eine stationäre Gasbewegung ohne Wärmeleitung und in Abwesenheit von 
Massenkräften. Dann existiert bekanntlich eine Grenzgeschwindigkeit v,, die konstant längs 
jeder Stromlinie ist und daher (W2) befriedigt. Da die Bewegung stationär ist, ist (W 4) für jede 
Ortsfunktion A von selbst erfüllt. 
ER Ferner betrachte ich nur die vor kurzem von R. Prim) eingeführten, sog. ‚„P-Gase“. 
Er Diese werden dadurch definiert, daß sie einer Zustandsgleichung von der Funktionalform 


N | o=P(p) H(n) (n= spez. Entropie)‘. ... ..... (35) 
Ei genügen. Sie enthalten offenbar die idealen Gase als Spezialfall. Nach R. Prim besteht das 


vollständige Gleichungssystem der stationären Bewegungen von P-Gasen unter den obigen 
Voraussetzungen aus der Kontinuitätsgleichung 


a. 


div (Pip) v0. a De (86), 
der Integrabilitätsbedingung für die Eulersche Bewegungsgleichung FERN. 
D.xX im, 
rot ee) N a 87 
und der Bernoullischen Gleichung | 
II(p) 
1 EN een TEE 38 
i II(p.) (38) 
Dabei ist //(p) definiert durch : 
Do) = UN Vs a DE Re : (39), 


?, ist der längs jeder Stromlinie Konstante Staudruck. Nach (36) kann man (W5) durch die 


_ Wahl » = P(p) befriedigen. Wegen (37) kann (W6) durch die Festsetzung 4 ua erfüllt 


II(p,) 


werden. Dann läßt sich infolge Gleichung (38) A auch schreiben A 270) ? 


Dayın folgt auf 
Grund unseres kinematischen Wirbelsatzes 


Der neue Wirbelsatz der Gasdynamik ! 


Gegeben sei eine stetige, stationäre, komplex-lamellare, massenkräftefreie Pr -Gas-Bewegung 
ohne Wärmeleitung. Dann ist 


% W, w. dlog II(p) | 
ER — a = — ee = c00nSt 2 2. nen 
"Up)Po) h dp | | 
längs jeder Stromlinie, d.h. ist eine angebbare Funktion des Druckes. Die Funktion h ist 
dabei nach (W1) definiert. e 
Dieses Ergebnis läßt sich mit dem schönen, von P. Nemen yiund R. Per 16) a I 


1947 entdeckten Wirbelsatz für Beltramische Gasbewegungen 102 x sn = N) yergleighen.”) un 
Im Spezialfall einer ebenen Bewegung wird aus (40) \ 


| % 1a 2ge 17 (oe = const ARE ei a A 
7 dp EN 
15) (1949, 1.] 16) [1949, 2, Theorem 6.] 


f 
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; Für ein ideales Gas gilt I/(p) > p * ‚ woraus unmittelbar mit Hilfe von (41) der erste Crocco- 
| Primsche Satz folgt. 
E . Im Spezialfall einer rotationssymmetrischen Bewegung wird aus (40) 
E- w, dlog II( 
; _. NE on Lea): 
| % dp 


woraus sofort der zweite Crocco-Primsche Satz folgt. 

Damit sind die Fragen 2 und 3 beantwortet. 

Da der neue Wirbelsatz der klassischen Hydrodynamik und der neue Wirbelsatz der Gas- 
dynamik Spezialfälle unseres allgemeinen kinematischen Wirbelsatzes sind, ist auch die Frage 4 

beantwortet. 

Es wurde schon bemerkt, daß der d’Alembertsche Wirbelsatz und der Svanbergsche 
Wirbelsatz unmittelbar aus den Helmholtzschen Wirbelsätzen folgen. Da wir jetzt eine enge 
Verwandtschaft zwischen dem d’Alembertschen Satz, dem Svanbergschen Satz und den 
beiden Crocco-Primschen Sätzen entdeckt haben, erhebt sich nun ganz natürlich die Frage, 
ob die letzteren (entsprechend der Verallgemeinerung der d’Alembertschen und Svanberg- 
schen Ergebnisse in den Helmholtzschen Sätzen) nicht auch verallgemeinert werden können. 
Leider ist die Antwort darauf negativ. In der klassischen Hydrodynamik ist, wie wir schon ge- 

E sehen haben, die Wahl A=1 eine zulässige; (18) wird dann die wohlbekannte ‚„‚Helmholtzsche‘ 
i Gleichung, aus der die Helmholtzschen Sätze ableitbar sind. In der Gasdynamik wird aber 
mit A= const die wichtige Bedingung (W6) nicht erfüllt; das letzte Glied rechts in (18) ver- 
schwindet nur, wenn A längs der Wirbellinien konstant bleibt. Daher ist Gleichung (18) im 
allgemeinen nicht die analoge Gleichung zur ‚„‚Helmholtzschen“, und es gelten deshalb keine 
"Sätze in der allgemeinen räumlichen .Gasdynamik, die den Helmholtzschen Sätzen direkt 
entsprechen. 

Ich vermute, daß noch andere wichtige Anwendungen des neuen kinematischen Wirbel- 
betrags-Satzes möglich sind, darunter vielleicht solche, die Bewegungen von zähen oder hetero- 
genen Flüssigkeiten betreffen. 

Die vorliegende Arbeit bildet einen Teil einer eingehenden Untersuchung der allgemeinen 
Kinematik von wirbelhaften Bewegungen 1”). 
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Lösung eines Minimumproblems der Tragflügeltheorie 
Von Karl Nickel in Tübingen 


Es wird die dritte Grundaufgabe der Prandilschen Theorie der tragenden Linie auf den Fall endlich 
vieler linearer Nebenbedingungen erweitert und gelöst. Drei Beispiele zeigen praktisch vorkommende Fälle 
mit solchen Nebenbedingungen. 


The third fundamental problem of Prandil’s theory of the supporting line is extended to the case of a finite 
number of linear accessory conditions, and the solution is given. Three examples are demonstraling Ihe 
practical occurence of such cases. 


Le troisiöme probleme fondamental de la theorie de la ligne de soutien de Prandtl est &tendu au cas d’un 
nombre fini de conditions secondaires lineaires; on en donne la solution. Trois examples montrent de tels 
cas d’oceurrence pratique. 

B 9Toit pabore oÖoömaerca u pemacrTca TpeTbA OCHOBHaA Banaya Teopun Upanaraıa 
AA CIyyaAd KOHEYHOTO'YHCHA JIMHEÜHBIX MONONHNTENBHEIX ycaoBnü. B Tpex Ipumepax 
HPHBONATCA BETpeyamımnecA Ha IIPAKTHKe CIyyam C TAKHMH AOUOAHHTEJIBHBEIMH YCHIOBHAMH. 


1. Aufgabenstellung 


In der Tragflügeltheorie I der „‚Vier Abhandlungen‘) gibt L. Prandtl als ‚Dritte Grund- 
aufgabe‘ das folgende Miniumproblem an (a. a.O., S. 28): 

„Gegeben ist der Gesamtauftrieb und die Flügelspannweite, außerdem o und V; gesucht ist 
diejenige Verteilung des Auftriebs über die Spannweite, bei der der Widerstand ein Minimum 
wird.“ 

Mit Hilfe der angegebenen Formeln (a.a.O., S.27) lautet die Kalten Fassung 
dieser Aufgabe in einem Koordinatensystem, in dem sich der Tragflügel von —1 bis +1erstreckt: 

Es sei 


1 karly) dy 3) 


Uhr an) dy ©—y 


-) 


und 7(—1)= T(+1)=0 [dabei ist T’(x) die örtliche Auftriebsdichte, w(x) die induzierte 
Abwindgeschwindigkeit an der Stelle x des Flügels]. Die Funktion I’ (@) möge im Intervall 
<—1, +1) so bestimmt werden, daß 


” | Ele) viane 


zu einem Minimum wird unter der Nebenbedingung 


EL) aa An ae 


A = Gesamtauftrieb beliebig vorgegeben). 


Diese Aufgabenstellung wurde von M. Munk?®) auf im Faume beliebig gerichteten und 
verteilten Auftrieb erweitert und allgemein gelöst. 

Man kann das obige. Minimumproblem auch noch auf ei: Art verallgemeinern, indem 
man statt der Nebenbedingung (3) andere — und evtl. auch mehrere — Nebenbedingungen vor- 
„schreibt. Drei Beispiele sollen das näher erläutern: 

a) Ein fliegender Tragflügel soll durch einen (kleinen) Querruderausschlag in eine (flache) 
Kurve gelegt werden. Welche Gestalt muß die Auftriebsverteilung besitzen, die durch den Quer- 
ruderausschlag zusätzlich zu der symmetrischen Verteilung hinzukommt ‚damit die Erhöhung des 
induzierten Widerstands möglichst klein bleibt? Kennt man diese Verteilung, so kann man sie 
durch passende Rudergestalt annähern, um so beim Querruderausschlag überflüssige Verluste 
zu vermeiden. Die Frage lautet hier also: Gesucht ist diejenige Afitriebsverteilung, die du 
induzierten Widerstand (2) zu einem Me: macht, wenn das Rollmoment 


n 


einen vorgeschriebenen Wert besitzen 6 — Die Verwendung von Gleichung a) ist in diesem 


Falle exakt richtignur im ersten Augenblick des Querruderausschlages, solange noch keine Dreh- _ 


und Rollbewegung vorhanden ist, da nur dann die abgehenden Wirbellinien geradlinig sind. 
1!) L. Prandtl und A. Betz: Vier Abhandlungen zur Hydrodynamik und Aerodynamik, Neudruck aus 


den Verhandlungen des III. Internationalen Mathematiker-Kongresses zu Heidelberg und BaB, an Nubeuiizi 


der Gesellschaft der Wissenschaften zu Göttingen. Göttingen 1927. |’ 
.) Unter dem Zeichen soll der Cauchysche Hauptwert des Integrals verstanden sein. 


1919. 


®) M.Munk: Tsoperiniehtiet :he Ban ben aus der Theorie a Eye ER ee Gattingen | | 
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Immerhin kann man für kleine Dreh- und Rollgeschwindigkeiten die Lösung des formulierten 
Problems noch als Näherungslösung verwenden. 

Eine andere Aufgabe, die-auf dasselbe Problem führt, und bei der diese Schwierigkeiten 
nicht auftreten, wurde mir freundlicherweise von Herrn Prof. Dr. Prandtl mitgeteilt. Es ist 
dies die Frage nach der Auftriebsverteilung geringsten Widerstandes für einen Tragflügel mit 
exzentrisch angebrachter Last, durch die ein Rollmoment (4) erzeugt wird. 

b) Die Lösung der von Prandtl angegebenen dritten Grundaufgabe lautet [vgl. !), S. 32] 
bei der obigen Formulierung: 5 
DE wer yı— 

’ 7 i 
d. h. der Auftrieb ist nach Art einer Halbellipse über die Spannweite verteilt (vgl. die ausgezogene 
Linie in Abb. 1). Ändert man nun /'(x) in der gestrichelt eingezeichneten Weise etwas ab, wo- 
bei A in (3) seinen Wert behalten soll, so wirdsich der induzierte Widerstand (in der Nähe seines 
Minimums) nur wenig ändern. Das Holmbiegemoment an der Flügelwurzel 


[r@) zdz, 


das man für symmetrische Auftriebsverteilungen (wenn von einem Faktor 1/2 abgesehen wird) 
durch 


erselzen kann, wird jedoch kleiner geworden sein. 

Für einen freifliegenden Tragflügel bedeutet das folgendes: Infolge der verringerten Bean- 
spruchung kann der Flügel gewichtlich leichter gebaut werden. Da aber im stationären Waag- 
rechtflug Auftrieb und Gewicht gleich groß sind, heißt das, daß auch der Gesamtauftrieb A in (3) 
kleiner wird. 


ix) 


Bild 1 Bild 2 


Ersetzt man I'(x) durch A T'(x), so ändern sich Auftrieb (3) und Holmbiegemoment (5) 
proportional zu A, während der induzierte Widerstand proportional zu A?ist. Man wird also mit 
einer Auftriebsverteilung, die an den Flügelenden etwas kleinere Werte als die elliptische Ver- 
teilung mit gleichem Auftrieb besitzt, etwas günstigere Widerstandsverhältnisse erwarten dürfen, 
als mit der elliptischen Auftriebsverteilung. — Allerdings muß bemerkt werden, daß die Ab- 
hängigkeit des Fluggewichts vom Holmbiegemoment bei den üblichen Tragflügelbauweisen sehr 
gering ist, so daß der oben beschriebene Effekt klein ist. 

Diese Überlegung führt auf die folgende Aufgabenstellung: Es soll der induzierte Wider- 
stand (2) zum Minimum gemacht werden unter den beiden Nebenbedingungen, daß Auftrieb (3) 
und Holmbiegemoment (5) vorgeschriebene Werte besitzen. 

 e) Betrachtet man einen in Flugrichtung gekrümmten Tragflügel (vgl. Bild 2), so wäre für 
‘diesen bei der Lösung der dritten Grundaufgabe .der Tragflügeltheorie 2) außer dem Gesamt- 
auftrieb auch noch das Längsmoment konstant zu halten. Man kann nun verlangen, daß die 
Erhaltung des Längsmoments ausschließlich durch Auftriebskräfte am Tragflügel erfolgen soll 
(Längsstabilisierung durch Tragflügelverwindung, wie etwa im Vogelflug). Sei etwa h(x) der 
Abstand der Druckmittelpunktslinie (in Bild 2 gestrichelt) von der x-Achse. Dann bedeutet 


diese Forderung, daß N 
DEE ZU 22. nel 
—1i 


_ (Längsmoment M vorgegeben) sein soll und die obige Aufgabe lautet: Gesucht wird diejenige 
Auftriebsverteilung /'(z), für die (2) zu einem Minimum wird unter den beiden Nebenbedingungen 
(3) und (6). — Es mag hier.angemerkt werden, daß in dem allgemeinen Falle eines beliebig ge- 


er war läßt sich (vgl. !), S. 25) das Bild der tragenden Linie auf diesen Fall nicht mehr anwenden, wie 
jedoch M. Munk (°), 8.21) gezeigt hat, genügt es bei Widerstandsbetrachtungen, nur das ebene Problem 
eines in Flugrichtung nicht gekrümmten Flügels zu betrachten. I'(z) ist hier also die Projektion des Auftricbs 
auf eine zur Flugrichtung senkrechte Ebene. 


werden: 


x h ; ee 17 | R EEE \ fi 
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stalteten Tragflügels die Druckmittelpunktslinie nicht die Linie des vorderen Neutralpunkts im 
ebenen Problem (7/4-Linie) sein wird. Die Aufgabe, zu einer beliebig vorgegebenen Tragfläche 
diese Linie zu bestimmen, bereitet große Schwierigkeiten, die m. W. bis heute noch nicht 
gelöst sind. Näherungsweise kann man jedoch die Druckmittelpunktslinie durch die 7'/4-Linie 
ersetzen, solange nur der Flügel nicht zu sehr von einem geraden Tragflügel abweicht. 


Es liegt nun nahe, auf Grund dieser Beispiele das oben formulierte Minimumproblem so zu 
erweitern, daß endlich viele, in /'‘(z) lineare Nebenbedingungen vorgeschrieben werden. Man 
sucht also (x) im Intervall <—1, +1) so zu bestimmen, daß 


1 ‘ 
SH 3 
zum Minimum wird unter den N Nebenbedingungen 
1 
J I(&) (2) do = 4A, a Be 
a ? 


[L—1)=T(+1)=0; h,(x) und A, beliebig vorgegeben, einzig für h,,(x)= 0 ist natürlich 
nur 4,0 zulässig). Für N=1, h,(x) = 1, A, = Aergibt sich daraus die von L. Prandtl 
formulierte Aufgabe als Sonderfall, ebenso lassen sich die Beispiele a), b) und c) durch analoge 
Spezialisierung gewinnen. 
- Umformung des Problems 
Zur mathematischen a des Problems mögen folgende er vorgenommen 


Man setzt 
HEICOSISE y= cost, 


4u@ fi —a=fl),  Tl)=Z(s) 


und wendet die Schreibweise %,(x) /1 — x? = h, 6) (n—1,..., N) an. Damit wird 5) aus (1), (2) 
(unter Vernachlässigung eines Faktors 1/4) und (N): 


1 "dzı) sin sdt 


° . 
Allen m. dt cos 1— cos s 8 e 
1 AO) OL (2a), 
[20 h„(s)ds= A, ee) 
Bei bekanntem f(s) ist (la) eine Fredholmsche se Re erster Art für Ihre 
Lösung ist bekanntlich ®) : Dt» 
20 _ 1 16 sin idt HE A 
| Mikes 9, EB 
woraus durch Integration 7% 
| | N na ” E a 2 
Zi): a. EN m 2 se a 
don) uch ni di en 
r sı D) I Fer 
folgt. Schreibt man zur Abkürzung £ ur: in 
i | ..s+t t Ä 
sin —,- BRUT Ta ar 
en log is == S (s, ) 5 " 3 j ei er vis 
sin —— BGE 79 
a: 2 ; 


Er 
5) Die Voraussetzungen, unter denen die vorkommenden Int Hals und Re A Beihen « e2 
und unter denen die verwendeten Umformungen (Vertauschung ı der DR ocratinnereitihfolge, 

Integration und Summation) erlaubt sind, werden in einer ausführlicheren Arbeit | 
unter dem Titel „„Lösung eines speziellen Minimumproblems“ in der Mathematischen 
- Außerdem wird darin noch eine zweite Klasse von Nebenbedingungen zu 

Zusatz bei der Korrektur: Ist inzwischen erschienen in der Math. Zeits 
6) Vgl. etwa K. Schröder: Über eine aan ax ‚erstei | 
berichte dr Preußischen ee der Wissenschaften XXX oe 
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und setzt (8) in (2a) und (7a) ein, sö erhält man 


/ J BONS EN TO (2b) 
_ und 
F 
I. | J I BS@ a de en 1...,. Nase. (2b). 
Da Burch (1a) und (8) die iosen Z(s) und f(s) einander eindeutig zugeordnet sind, kann 
man auch /(s) als gesucht ansehen. Damit lautet die zu lösende Aufgabe: Es soll f(s) im Intervall 
<0, ©) so bestimmt werden, daß (2b) unter den Nebenbedingungen (7b) zum Minimum wird. 
Zur weiteren Vereinfachung der Schreibweise werde noch folgende Abkürzung eingeführt: 
G )=/ IRIORIOEIG RR re (9). 
| Mit der für 0 < st In gar gültigen Darstellung’) | 
2 - 
Slate * 2 sin nssin nt 
| ergibt sich aus (9): 
w (9) = In (s) g(t) S(s, 1) ds dt 
} x 9 oo 
yes an u /(s)sin neu) g(t) sin nid). 
| oder = 
(fN)= Pr a IE SEEN BE TE, (10) 
i mit den Fouriersinuskoeffizienten 
u IR 
a, -2 | f(s)sinnsds, | g)sinnid (n=1,2,...). 
| 0 0 j 
| Nach (9) und (10), und wegen S(s, t) = $(t, s) besitzt (f,g) die folgenden Eigenschaften: 
3 ® (af, g) = a V g) (a= reelle Zahl) 
| = HEN EN + A 
dd EN>O für Fee) =:0 
\ Man hat also die Aufgabe, eine Funktion f(s) im Intervall <0, x) so zu finden, daß 
: (f, f)= Min. | 


wird unter den N Nebenbedingungen  - B 


E*. (, h,)= 4, Melt, Ni lan... (12). 
Dabei ist (/, g) durch (9) erklärt und hat die Eigenschaften (11). 


3. Lösung des Problems 
Wenn die Funktionen 7, (s) linear abhängig sind, so gibt es unter den Nebenbedingungen (12) 
entweder überzählige, oder sie sind unter sich widerspruchsvoll, lassen sich also überhaupt nicht 
erfüllen. Man darf also, ohne die Allgemeinheit einzuschränken, die Funktionen h,(s) als linear 
unabhängig voraussetzen. 
Borsuasesgung: Die Funktionen A,(s) seien linear unabhängig, die reellen Zahlen 
> An ag 1, ...., N) seien beliebig gewählt. 
Es; npräng- Dann gibt es genau eine Lösung f’(s) des obigen Minimumproblems und 
# % diese Lösung besitzt die Gestalt 


RR | FO=L a6) 
it einde ei bestimmten. reellen Zahlen a,. 
etwa G.Hamel, Integralgleichungen 8. 120. Berlin 1937 bei J. Springer; sowie K, 


Ermittlung einer Reihendarstellun Ne. er In r in elliptischen Koordinaten. Z,. angew. 
“4 ur 30 (soo) Seite 186 Formel (16). 


a ap 1 Br N U A Ee 
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Beweis: Nach dem Orthonormalisierungsverfahren vonE. Schmidt kann man Funktionen 
= & Onnfin (s) (m FÜR EN) 


M= 


so finden daß 1 fü mn 


0 O fü m#n 
ist. Als Umkehrung dieser Gleichungen ergibt sich eine Darstellung 


h„(s) === 53 dm Anm (8) (d,, + 0; Hz 1, a N). 
m=1 


NER H,,) — On 


Man setzt nun 


5 
fo) =2 „H,(s) 
n=1 
mit vorläufig noch unbestimmten el en c„ und erhält mit den Rechenregeln (11) 


#3 > (BER) > €, 2 ER PR H,,) 


n—=1 


N r 
pH On I Umr® = Dt für. m 1 Bee 


n=1l m=1 n=1 x 
Sollnun I c,d,,=4, (r=1,...., N) sein, so sind dadurch die c, rekurrent eindeutig be- 
nl 


stimmt. Mit den so gewählten reellen Zahlen c, erfüllt dann f’ die Nebenbedingungen (12), weiter 
erhält f’(s) die behauptete Be 


ih (s) = EI ©) Cm An (5) = Se 2 Cm = Cnm h, (= > h, (3) - 5 m =: I Ay h (8). 


Es bleibt Boch, zu zeigen, daß Fe Ri (f, f) zum Minen? ERNS Nun läßt sich eine be- 
liebige Funktion f(s) immer in der Gestalt 


FF) HKE) mämlich mit KA) fl) Fl) 
schreiben. Damit ist I) 
(7 h,) Fr (er k, h,) = ee h,) m (k, h,) =4A,+ (k, h,) (n = l, ..., N). 
f(s) erfüllt also die Nebenbedingungen (12) genau dann, wenn 
(RO isWwlüreeNn EIER ENR 
Für die Gesamtheit der Funktionen f(s), die die Nebenbedingungen (12) erfüllen, gilt also 


nach den Regeln (11): 
AN a EEE) Tal EM 
a + 1) =(f P) 


wobei das Gleichheitszeichen nur für k(s) =0 eintritt. Die Eindeutigkeit der Koeffizienten a, 
folgt mit (11d) aus der linearen Unabhängigkeit der Funktionen A,(s). Damit ist die Be- 
hauptung "bewiesen. 
Macht man die oben eingeführten Transformationen wieder rückgängig, so findet man 
folgendes Ergebnis: 
4. Ergebnis 


Es seien die Funktionen h, (x) und die reellen Zahlen A, (n=1,..., N) beliebig gewählt. 
Die Funktion Z'(z) sollin —1<x=<-+1 so Br werden, daß Ts) —= T(+1)=0ist 
und daß unter der Setzung 


a CAr(y) dy 
| w(a) = na ae (1) 
3 A 1 
der Ausdruck J I(o) wie) dal. nos ae (2). 
SEN ) 4 


zu einem Minimum wird unter den DER neingElE En ; 4 . 
Jo) 1.(a) de= 4, a rn Wim Are een 


Wenn es überhaupt eine Funktion I (x) gibt, die den N ee (7) Benieh dann 
hat dieses er Ba eine Lösung, und diese besitzt die Gestalt: 


I(a)= = a, «| h(y) log een N I R 


% Ba 
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mit reellen Koeffizienten @,. Sind die Funktionen h,(x) linear unabhängig, so hat das Minimum- 
problem immer eine Lösung und die Koeffizienten «a, sind eindeutig bestimmt. 


5. Anwendung 


In der folgenden Tabelle stehen zu einigen einfachen Funktionen Ah(x) die zugehörigen 
1 


— a) — 2) 
Funktionen | h(y) De a El N dy. 
Br HTERd, 
+1 2 
Kia) | ring HH ee A 
ie 2 —y 
1 ay1 — a 
x Sala 
5 eaınjioa 
23 ; 2 +1)j1— a2 
Er 9 een 
a] r «| 
Han Bar 

el x log a 2 -/i—2 

2 1 1— in 
2 |e| > (@108 En, Et ale) 

1 1— 2 — 
ze] 1 (Barton FE a nyize) 
Ben zaz $: 1 ma + Al — ©) — 
ne (2, — x) log er, +Y1—a2arccos z, 


Damit findet man für die ersten beiden der eingangs angegebenen Beispiele die Lösungen: 


Beispiel Nebenbedingungen Lösung | 
1. U NERLIE 
a) [ Te) zd2= R IK zyl— 2 
1 It 
9 LER ER 
[r@ar-4| r@=3(? _a)yr—a 
1 
b) 
[f@) \elde=H +3 (38 — ra = 
= 


In dem Bild 3 ist diese Lösung des Beispiels a) aufgetragen. Bild 4 zeigt einige Kurven 
aus der Schar der Lösungen von Beispiel b) für A = constg H variabel (unter diesen ist für 
3rH—=4A4 die Kia hellipes als Sonderfall enthalten). 


Too 


Bild 3 zu Beispiel a) Bild 4 zu Beispiel b) 
Eingegangen am 28. Januar 1950. 
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Finige exakte Lösungen für die Temperaturverteilung 


in einer laminaren Strömung 
Von H. Schlichting in Braunschweig 
Für die ebene Kanalsirömung und für die Strömung zwischen zwei konzenirischen rotierenden Zylindern 
werden ewakte Lösungen für die Temperaturverteilung mit Berücksichtigung der Reibungswärme angegeben. 


For ihe plane canal flow and the flow beiween two concentric rotaling eylinders, the distribution of tem- 
peralure is given wilh regard to the friectional heat. 


Pour l’£coulement plan dans un canal et pour V’&coulement entre deuz cylindres nen rolatives 
on donne la distribution ewacte de la temp£rature en consid£ration de la chaleur de friction. 

Jaercs Toynoe pemernme sanaym 0 pacnpeneneHnn TeMNepaTypH C YYeToM TEIAOTH 
TpeHunA B CIyyae INIOCKOTO TeYCHNA YKHNKOCTU B KAHAJle MHIH TEeYeHNHA 3EHNKOCTU MEiKAy 
AByMA KOHHeHTPHYHBIMH BPAILOIHMUCH IMIHHNPAMH. 


I. Einleitung 


Für die Fragen des konvektiven Wärmeüberganges zwischen einem festen Körper und einer 
strömenden Flüssigkeit (oder einem Gas) spielt die Temperaturverteilung in der Umgebung des 
umströmten Körpers eine grundlegende Rolle. Die genaue Ermittlung dieser Temperaturverteilung 
ist im wesentlichen ein Strömungsproblem. Bei einigermaßen großen Strömungsgeschwindig- 
keiten hat das Temperaturfeld ‚„‚Grenzschichtcharakter“, d.h. es vollzieht sich der Temperatur- 
ausgleich zwischen der Flüssigkeit und dem umströmten Körper in einer sehr dünnen Schicht, 
der sog. ‚„ Temperaturgrenzschicht“, deren Dicke sehr klein ist gegen die Abmessungen des Körpers. 
Bei Gasen ist die Dicke der Temperaturgrenzschicht von etwa gleicher Größe wie die Dicke der 
Strömungsgrenzschicht, während bei Flüssigkeiten die Dicke der Temperaturgrenzschicht u. U. 
erheblich geringer ist als diejenige der Strömungsgrenzschicht. 

Für die Temperaturverteilung in der Umgebung eines geheizten umströmten Körpers sind 
in den letzten Jahren verschiedene Näherungslösungen auf der Grundlage der Grenzschicht- 
theorie angegeben worden!-®). Die Zahl der bekannten exakten Lösungen ist jedoch sehr 
gering. Sie beschränken sich im wesentlichen auf die laminare Kanal- und Rohrströmung?°), und 
hierbei auf solche Fälle, wo die in der Strömung durch Reibung erzeugte Wärme vernach- 
lässigt wird, was nur für mäßige Strömungsgeschwindigkeiten zulässig ist. Im folgenden 
sollen einige exakte Lösungen für die Temperaturverteilung in laminarer Strömung mit Be- 
rücksichtigung der durch Reibung erzeugten Wärme angegeben werden. Diese Lösungen haben 
einmal eine gewisse grundsätzliche Bedeutung, da sie exakte Lösungen der Strömungsgleichungen 
und der Gleichung für die Temperaturverteilung darstellen. Andererseits haben einige davon 
auch eine unmittelbare praktische Bedeutung für die Theorie der hydrodynamischen Lagerreibung. 
Trotz ihrer sehr großen Einfachheit scheinen die meisten dieser Lösungen für die Temperatur- 
verteilung bisher noch nicht bekannt zu sein. 


U. Die Grundgleiehungen 
Die Strömungsgleichungen (Navier-Stokessche Gleichungen) und die Gleichungen für die 
Temperaturverteilung (Energiegleichung) lauten für ein inkompressibles Medium mit tempera- 
turunabhängigen Stoffbeiwerten für den ebenen Fall in kartesischen Koordinaten: 


= ee, 
rt) @ 
olu 2” te). N ee 
(u 22 4037) = tal rRB 


Körper. VDI-Forsch.-Heft 416 (1942). 


®) E. Eckert und 0. Drewitz: Der Wärmeübergang an eine mit großer Geschwindigkeit TAngESnBenE # 


strömten Platte. Forschung a. d. Gebiet d. Ingenieurwesens, Bd. 11 (1940), 8. 116. 
®) E. Eckert: Wärme- und Stoffaustausch. Berlin: Springer 1949. 
*) J. Ginzel: Ein Pohlhausen-Verfahren zur Berechnung laminarer kompressibler ie an 
einer geheizten Wand, Z.angew. Math. Mech. 29 (1949), S. 321. 
®) M. Jacobs: Allgemeine Grundlagen der Wärmeübertragung. Der OhehnitrIngeisesl| L Teiprig us 


E. Eckert: Die Berechnung des Wärmeüberganges in der laminaren Grenzschicht umströmter > 
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Dabei bedeuten «, v die Geschwindigkeitskomponenten, 7 die Temperatur, p den Druck, o die 
Dichte, v[kgs/m?] die Zähigkeit, 9 die Erdschwere, c,[m/grd] die spezifische Wärme bezogen 
auf die Gewichtseinheit, A[kg/s grd] die Wärmeleitfähigkeit und ® die Dissipationsfunktion, 


welche gegeben ist durch 
9u\? (@)\ ( E) 
TH re er 2 RE edle 
H E ar ara er 6): 


In der Gl. (4) für die Temperaturverteilung bedeutet die linke Seite die Wärmezufuhr pro Volumen- 
einheit und Zeiteinheit durch Konvektion, die ersten beiden Glieder der rechten Seite die Wärme- 
zufuhr durch Leitung und das letzte Glied die durch Reibung erzeugte Wärme. Der Wärme- 
übergang zwischen der Flüssigkeit und den Wänden ist gegeben durch den Wärmestrom q 
(= Wärmemenge pro Zeiteinheit und Flächeneinheit). Es gilt 


Be a. LEN. (BE 


wo 2T/du den Temperaturgradienten an der Wand bedeutet. Dabei ist q> 0 für einen Wärme- 
übergang von Flüssigkeit auf Wand, wenn die Normale n positiv gerechnet wird in Richtung 
Wand — Flüssigkeit. 
- III. Couette-Strömung | 
Eine erste exakte Lösung für das obige Gleichungssystem läßt sich angeben für die 
| Couetteströmung zwischen zwei parallelen ebenen Wänden (Wandabstand 3), von denen die 
eine in Ruhe ist und die andere sich mit der konstanten Geschwindigkeit U bewegt. Die Rand- 
bedingungen sind somit für die Geschwindigkeit 


y ER Al und WERIYU—il; 0: 


Die Randbedingungen für die Temperatur sind, daß entweder die Temperatur beider Wände 
vorgeschrieben ist (Wärmeübergang an beiden Wänden) 


Re Tan N N za 


oder daß nur an einer Wand die Temperatur vorgeschrieben ist, während an der anderen Wand 
verlangt wird, daß kein Wärmeübergang stattfindet. In diesem Fall findet dann eine Auf- 
heizung der letzteren durch die Reibungswärme statt (Thermometerproblem): 


oT 
Die Lösung der Strömungsgleichungen, welche die Randbedingungen erfüllt, lautet: 
u-U4; ZUR Dec RDUSUr ae ar (Oh 
Eine von x unabhängige Lösung für die Temperaturverteilung erhält man dann nach (4) aus 
2 ee Me 
Ä dp #\ay 
oder mit GI. (8): 
@T (4 09 
a 2 ak BR Nr 
Die Lösung von Gl. (9), welche den Randbedingungen (7a) genügt, lautet 
i BY ER, 9 ea 4) 
T(y)— Tı=(T; 7), + 31 gr 1 n D . . . . . . . (10). 


Hiernach überlagert sich der linearen Temperaturverteilung in der ruhenden Flüssigkeit 
eine parabolische Verteilung, die von der Reibungswärme herrührt. Die Temperaturverteilung 


S } $ U° ER 
gi für verschiedene Werte der Dimensionslosen en AL 7 Rn) ist in Bild 1 dargestellt. 


Wärmeübergang: Für den Fall, daß beide Wände wärmedurchlässig sind, erhält man 
‚aus (6) und (10) für den Wärmestrom | 


| | Ip 0 

ander unteren Wand (u=0): ee h +5, a HAFEN: 
br FE, ur ! . 4 T No T 12 

= ander oberen Wand w=D): = +, KR ib). 
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Der gesamte Wärmestrom an beiden Wänden, g=g, 4-9, der im Falle der Ruhe (oder bei 


Vernachlässigung der Reibungswärme) verschwindet, ist im allgemeinen Fall: 
uU? 
Ta 


also unabhängig von der Tem- 
peraturdifferenz der Wände. Be- 
merkenswert ist, daß bei vor- 
gegebener Temperaturdifferenz der 
beiden Wände, 7,— T,>0, ein 
Wärmeübergang von der wärmeren 
Wand auf die Flüssigkeit (9,< 0) 
nur eintritt, solange die Geschwin- 
digkeit der wärmeren Wand unter- 
‚halb einer bestimmten Grenze 
bleibt. Man erhält Kühlung der 
wärmeren Wand (9, < 0) bzw. Auf- 
heizungder wärmerenWand (9 > 0), 
«U 


L 
Bild 1. Temperaturverteilung bei der Couetteströmung mit Berücksichtigung Je nachdem ob S (7, SE, IHR 
der Reibungswärme, © = „U?/A (T,—T,). sul 21 


Im Falle, daß beide Wände gleiche Temperatur haben (7, = T,), erhält man infolge der 
Reibungswärme eine parabolische Temperaturverteilung über die Spaltbreite mit der maximalen 
Temperaturerhöhung in der Mitte vom Betrage 

s | 
nn N N 
Bemerkenswert ist, daß diese Temperaturerhöhung infolge der Reibungswärme unabhängig ist 
von der Spaltbreite A®). 
Für die Randbedingungen (7b) (Thermometerproblem) erhält man aus (9) die Lösung 


SE a fi ee) 
(y) A — Ta): 9, 66 Ve (13) 
und hieraus für die Aufheiztemperatur der unteren Wand 
u U2 1 
RE En ie (14). 


Diese Aufheizung (Eigentemperatur) beträgt also das Vierfache der maximalen Temperatur- 
erhöhung in der Spaltmitte für den Fall, wo beide Wände auf gleicher Temperatur gehalten werden. 

Die Temperaturverteilung nach Gl]. (10) und Bild 1 ist von Bedeutung bei der Strömung 
im Schmierspalt zwischen Zapfen und Lager. Diese ist bei den kleinen Abmessungen des Schmier- 
spaltes und der großen Zähigkeit des Schmieröles lJaminar. Die Aufheizung des Schmieröles durch 
die»Reibungswärme ist hier auch schon bei mäßigen Geschwindigkeiten beträchtlich, wie das 
folgende Zahlenbeispiel lehrt: 


Ölzähigkeit (bei etwa 30°C): u = 0,040 kg s/m?, 
Wärmeleitzahl: = 0,02 kg/s grd. 


Hiermit ergibt sich aus Gl. (12) für U=5m/s eine Temperaturerhöhung von T,— T,—= 6°, 


und für U= 10 m/s, T,„— T,= 25°. Diese Temperaturerhöhung des Öles ist so beträchtlich, 
daß sie bei der starken Temperaturabhängigkeit des Zähigkeitsbeiwertes für die Tragfähigkeit 
des Lagers wesentlich wird. 


IV, Kanalströmung 
Eine zweite recht einfache exakte Lösung für die Temperaturverteilung erhält man aus 


Gl. (1) bis (4) für die Poiseuillesche Kanalströmung. Mit den Bezeichnungen nach Bild 2 
ist die Geschwindigkeitsverteilung i 


uly) = u, (1 _) ’ 


6) Dieser Sonderfall ist bereits von G. Vogelpohl [VDI-Forsch.-Heft 425 (1949)] behandelt worden. N 
Darüber hinaus gibt Vogelpohl die Lösung für die Temperaturverteilung im „Einlauf“, d.h. für diejenige 
Spaltströmung, bei der im Eintrittsquerschnitt eine über die Spaltweite konstante Temperaturverteilung vor- y 


handen ist, die weiter stromabwärts allmählich in diese parabolische Verteilung übergeht. 
\ r B 


= » 


4 
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Dabei ist die Maximalgeschwindigkeit in der Kanalmitte «,, gegeben durch den Druckgradienten 


t 2 
längs der Kanalachse, «, = es) - Hiermit ergibt sich für die Temperaturverteilung 
F aus (4) und (5) eine nur von y abhängige Lösung aus 
4 @T 4 uu,, 
de zu gr y? Ca Su Fr (15). 


Wärmeübergang an beiden Wänden: Schreibt man für die Temperatur als Rand- 
bedingung vor T=T, bei y=—hund T=T, bei y=-+-h, so ergibt sich aus (15): 


1 1 A 
WESER)... (0). 


Das letzte Glied stellt die Temperaturverteilung infolge der Reibungswärme dar. 
Für den Wärmestrom erhält man aus (6) und (16): an der unteren Wand (T= T,}): 
: ; BAR = 2, 4 Kum 
2 Tan 


4 
und an der oberen Wand (T = T,): 


p=—4 


] Der gesamte Wärmestrom (an bei- 
den Wänden) ist 


Die Umkehr des Wärmeüberganges 
an der wärmeren Wand (% = 0) 
tritt in diesem Fall ein für 
sW 
3 u 7 = Mertae Wh 1° 
Für den Sonderfall, daß beide 
Wände auf gleicher Temperatur 
gehalten werden, 7,=T;=T,, er- 
gibt sich eine Temperaturvertei- 
lung nach einer einfachen Parabel 4. 
Grades, die in Bild 2 als Kurve a 


d . v c Bild 2. Temperaturverteilung bei der Poiseuilleschen Kanalströmung 
argestellt ist. Die maximale Tem- @ beide Wände haben gleiche Temperatur; Wärmeübergang an beiden Wänden, 
peraturerhöhung in der Kanalmitte db untere Wand wärmeundurchlässig, Wärmeübergang nur an der oberen Wand 


ist in diesem Fall (Thermometerproblem) 


- Auch hier ist die maximale Temperaturerhöhung unabhängig von der Kanalbreite.”) 
Thermometerproblem: Wird Wärmeübergang nur an einer Wand zugelassen, dagegen 


_  _  dieandere alswärmeundurchlässig angesehen, so lauten die Randbedingungen für die Temperatur: 
e: i _ T=T, bei y=+h und s; —0 bei y=—.h. Die Lösung für diesen Fall lautet: 

Bo | Be REN 7 
ER \ i Iy)—-T=3%5Wn\d 47 BR 0 OU (17). 


N N ET \ 7 u 
Fall beträgt also die Aufheizung der wärmeundurchlässigen Wand das Achtfache der 
ıalen Temperaturerhöhung in der Kanalmitte für den Fall, daß beide Wände auf gleicher 
ur gehalten werden. | 

ei ‚der Korrektur (Sept. 1950): Eine Erweiterung der Lösung für die Poiseui lle sche 
ine mit der Temperatur veränderliche Zähigkeit ist inzwischen von H. Hausenblas, 
151—166, angegeben worden. 
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V. Konzentrisehe rotierende Zylinder 


Eine dritte exakte Lösung läßt sich angeben für die Strömung zwischen zwei konzen- 
trischen rotierenden Zylindern, die mit verschiedener Winkelgeschwindigkeit umlaufen. 
Die Gleichungen mögen in ebenen Polarkoordinaten r, g angegeben werden (u Geschwindig- 
keitskomponente in der Umfangsrichtung). Da in diesem Falle alle Größen von @ unabhängig 
sind, reduzieren sich die Strömungsgleichungen auf 


u dp, \ deu A )= h % & 
en =; tz) Re 


Die Gleichung für die Temperaturverteilung 7(r) lautet, da die konvektiven Glieder wieder 
fortfallen, 
amt 21 EN 
R 0=la + or + uDa PEN (20) 


mit der Dissipationsfunktion 


Seien r,, r, die Radien des inneren bzw. äußeren Zylinders und ®,, &, die entsprechenden Winkel- 
geschwindigkeiten, so lautet die Lösung für die Geschwindigkeitsverteilung 


j 
E 
R 
} 
NE 


u(r) = pp: a ER, U a et cn Se. (22) 
mit | 
ana, 2 nenn 
r, —tr 7, — 
2 1 2 "8 


Aus Gl. (20) ergibt sich als Lösung für die Temperaturverteilung ohne Berücksichtigungder 
Reibungswärme mit den Randbedingungen T=T, für r=r, und T=T, für r=r;: 
5 h I 
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Bild 3. Temperaturverteilung in der Strömung zwischen 
zwei konzentrischen rotierenden Zylindern, ohne Rei- 


bungswärme. k = in 
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Mit Berücksichtigung der Reibungswärme ist nach‘ (20), (21), (22) die Gleichung 


aT u u b? R 
| Ele een mHn 222.0: 
zu lösen. Für die Randbedingungen 7T= T, bei r=r, und r = r,.ergibt sich 
4 In &) 
EL ENG ee A al 
T—T, FL Le A n(e)" oA 
2 2) 
Mit U, — (w, — ®,)r, sowie k— "2 und 7 — ”. Jäßt sich Gl. (26) in der dimensionslosen Form 
schreiben Mi 7] 
N ER 1 iR ‚en 12 =) 
; Fr ale ln a > a). 
u: 
Die hiernach berechnete Temperaturverteilung ist in Bild 4 für verschiedene k = 5 dar- 
gestellt. Im Grenzfall (r n ie FA 0, d.i. k—1, erhält man natürlich wieder die parabolische 
) Temperaturverteilung der Couetteströmung nach Gl]. (10), wenn man dort 7,= T, setzt. Die 


allgemeine Lösung mit Berücksichtigung der Reibungswärme, aber mit verschiedenen Wand- 
temperaturen am inneren und äußeren Zylinder, erhält man durch Überlagerung der Lösungen (24) 
und (26 

et Weise wie bei der Couetteströmung und bei der ebenen Kanalströmung läßt 
sich auch hier leicht die Lösung angeben für den Fall, daß ein Wärmeübergang nur auf dem 
einen der beiden Zylinder stattfindet. Auch lassen sich in jedem Fall die Ausdrücke für den 
Wärmestrom auf Grund der angegebenen Formeln sofort aufschreiben. 

Eingegangen am 21. Februar 1950. 


Die mathematischen Grundlagen zu einer Meßmethode 
des Schubmoduls zäher Flüssigkeiten 
Von Erich Pieruschka in Kassel!) 


Der Schubmodul zäher Flüssigkeiten kann mittels eines Schwingmeßgeräts bestimmt werden. Es werden 
hierzu eine Phasenverschiebungs- und eine Amplitudenmethode entwickelt. Weiter werden Versuche aus- 
gewertet, die mit diesem Gerät an der Technischen Hochschule Danzig 1942 gemacht wurden. Dabei zeigt 
sich, daß es Flüssigkeiten gibt, die dem Newton schen Stoffgeseiz (G — ©) nicht gehorchen. 


To measure ihe modulus of shearing of viscous liquids, an oscillation instrument can be used. A phase 
# hi difference method and an amplitude method are developed for that purpose. Besides, the author gives 
Br aninterpretation of euperiments, made in 1942 by an instrument of that sort at ihe Technische Hochschule 
« Danzig. It appears, that some liquids do not obey N ewton's law of material (G — »). 


Le module de eisaillement d’un fluide visquex peut Eire determine & Vaide d’un appareil de mesure oscil- 
latoire. A cette fin une methode d’amplitude ei une methode de decalage de phase sont developp£es. En outre 
des esp£riences sont interprötees, faites avec cet appareil en 1942 @ la Technische Hochschule de Danzig. 

Il en resultait V’existence de fluides desobeissantes a la loi de substance de Newton (G> ). 


Ba, Monyas caBura BASKUX KHNKoCTeH MOKET ÖBITB ompezenen mp moMmoimu ImpaÖopa, ocHo- 
Pas BAHHOTO Ha IPHHuNNe BuÖpanum. Ü 9TOU MeNBIO aBTOPOM paspadartsıBaeTch MeTon cABura epas 
25: MH METON aMIMIHTyA. BareMm HMCcIeNyWTcA Pesy.IBTaTkl WCIBITAHHÜ, IPOHSBENCHHHX TIPu 
_  ToMoIm aroro ıpudopa » Bsremem Texnuuyeckom Vunsnmme B Nannre B 1942 rony. Ilpu orom 
Ka BBLACHAETCH, 9TO HEKOTOPHE HHAKOCTH He HORYET Bakony Marernm Hsiroma (@ > ©). 


I. Vorbereitung 


A r 1. Problemstellung 
in = | el ap Faksshem der Technischen Hochschule Danzig hat Herr Prof. Dr.-Ing. 
\ Hans] Fromm Messungen zur Bestimmung des Schubmoduls zäher Flüssigkeiten veranlaßt, die 
* lerrn wi ‚LELBUNG ausgeführt und ausgewertet wurden (beendet wurde diese Arbeit 1942). 


so sind N sschepretäkofle und Apparatur nicht mehr da. Nur einzelne Daten und Ergeb- 
| ten rekonstruiert werden, über die im wesentlichen H. Fromm?) in seiner Arbeit 


a tion an der Universität Mainz 1949. Berichter: Prof. Dr.-Ing. H. Fromm, Mitberichter: 
.» K. Bechert. T 
KR, 2. Aigen Math, A 28 (1948), 8. 43—54. 
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Zylinder soll sich nach (2) bewegen. 
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„Laminare Strömung Newtonscher ‘und Maxwellscher Flüssigkeiten“ berichtet. Hier, wird i 
auch auf Unstimmigkeiten hingewiesen, die sich damals nicht erklären ließen und ihre Ursache 
im Maxwellschen Stofigesetz 


TR 
en OL a 


das damals allein vorlag und eventuell auch in der Annäherung mathematischer Ansätze haben 
konnten. Nun hat Fromm in der oben erwähnten Arbeit das Maxwellsche Stoffgesetz durch 
eine Tensordarstellung, bei der auch die Drehung neben der Schiebung berücksichtigt wird, 
wesentlich verfeinert. Deshalb nahm ich die Anregung von Herrn Prof. Fromm auf, unter 
Verwendung des neuen Stoffgesetzes die genauen mathematischen Grundlagen zu der Danziger 
Versuchsanordnung zu schaffen und die Ergebnisse der Hartungschen Auswertung zu 
überprüfen. ; 


2. Die Differentialgleichungen des Problems 


Die Danziger Apparatur bestand im wesentlichen aus zwei konzentrischen Zylindern, die 
durch einen schmalen Ringspalt, der zur Aufnahme der Versuchsflüssigkeit diente, getrennt 
x waren. Der eine Zylinder wurde in Schwingungen konstanter 
Amplitude versetzt, die durch die Elastizität und Zähigkeit der 
Versuchsflüssigkeit auf den zweiten freischwingenden Schwing- 

ö zylinder übertragen wurden?). q 


y Bild 1 zeigt die Wahl des Koordinatensystems im Ring- 
spalt von der Breite d. Da diese gegenüber dem Durchmesser 
des Zylinders klein ist und ebenso die Amplitude des antreiben- 
den Zylinders nicht zu groß ist, kann man mit kartesischen 
statt Zylinder-Koordinaten rechnen. Der Weg des antreibenden Zylinders kann durch 


W=r-arcogat en. a. er u) 


Bild 1 


dargestellt werden. 

Bezeichnet man die Geschwindigkeit in y-Richtung mit v und die Relaxationszeit n/@ 
mit 7, so erhält man aus den Gl]. (7) der erwähnten Frommschen Arbeit für den Beharrungs- 
zustand unseres Problems 


:®) 
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Hier ist o=0,=—o, und 7=[7r,,. - 
Die dynamische Grundgleichung für die Schub- und Trögheiskntt in y-Richtung ah 
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wenn oe die Dichte der Flüssigkeit sein soll. _ RR Mr: 
Die Randbedingungen für den Beharrungszustand lauten: | A 


(0, u en 


\ v(ö, )= = silr Er elle * nt “ es SEE a REN: % 


wist die Flächenbelegung des freischwingenden Zylinders mit reduzierter Masse. De: antreibende BER 


Das Randwertproblem partieller Differentialgleichungen ) bis ©) gibt den 
vorgang vollständig wieder. Da die Dikterentialgisjetungen, zum Tei ie 
eine schrittweise angenäherte Lösung möglich. ” ie ale 

Bis zu einer gewissen Größe der Frequenz des RER; Zylin r: 
vernachlässigt werden, dann braucht man (4) nicht und die Fro mmscher ‚Glei 
Beh in den vereinfachten Ansatz (1) über, wie es min sein m 


: Eine genaue Beschhaikntn der Abpdralg und des Arbeiens 
Messung des Schubmoduls zäher Flünsigköiten. A 1951. 
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gewöhnliche lineare Differentialgleichung mit konstanten Koeffizienten, deren Lösung sehr ein- 
fach ist. Die erste wohl wesentliche Verbesserung der Hartungschen Lösung erhält man also, 
wenn man zwar o/@ vernachlässigt, jedoch o nicht Null setzt. Dieses Problem wird zuerst 
$ behandelt. 


II. Durehführung für das linearisierte System 
3. Lösung der Differentialgleichungen 


BA . 
# Bei Vernachlässigung von o/@ neben 1 erhält man aus (3) bis (7) zur Berechnung von 7 
Ya CRe lineare Randwertproblem: 
Dae '.@.08% 
rn er Re ec) 
a) 
(0, den E (9), 
? | | nn va RE Eye ale (0): 
; Dieses bedeutet Rechnung mit dem vereinfachten Maxwellschen Ansatz (1). 
a Dieses kann mittels des Ansatzes erfüllt werden: 
RN ; "r TER Smart ale) ©cosatn a (ddl) 
Er Setzt man a in 1 (8) bis (10) ein, so geht dieses Randwertproblem über in 
E: ° RIES 2 "tg otgr g(a@)=0 BT CR ARE RA) 
e ERERGE ORT: (8) — - ft LE REN ET ER 
mit den Randbedingungen: 
E Er =; JO) =, FO RL, THLLR 
Be; BR ER aa;  F=0 N RB (15). 
E ee ie Randwertproblem kann man'nach der Methode der linearen homogenen gewöhn- 
 liehen Differentialgleichungen lösen. Man erhält so für f(x) und g(2): 
Br: - Rae Je) = Kali (4} sin + A* cos va) esı®(A*sinvc+ Al cos vz) net (16), 
Br ‚ge = 2 ei Aayer — 4fcos va) + e 4° (— Afsinvc+ Afcosve).. . . (17), 
a wor N nt 
#0 DEE BEN, u H,=%9 und nn a (18) 
Rt 9 Et, O7 Pl „= ER und #®—=yl1l+9—)9 ER ENERE IT LONNS, (19) 
SERIE RT EHER: h 16 EASA,. i 
und Ro) £, 2 in EIN are Lana 
€ er ir u ee lg nn, ap DR ER ER N ar (20): 


Dart ie. ie Konstanten 4 En bis a, ergeben sich aus dem Gleichungssystem: 
za de a4 (di Me =, At. I ERARERERN 


\ 


ER A MR re (22), 

u sin vö+ a + I cos vö] = | (23) 
4) sin a A ) cos vö] 
Arysin »ö-r (Ag: + 42)cos völ + | (4). 


Br — AR) cos vö] = "a0? es? | 


\ jon en Vom besonderen Interesse ist hier nur 


86 Pieruschka, Die Aal einer Meßmethode des Nahe 8: 31 ne 2 ir 5: 
a VE EEE SEE EN I ee a er N, “ 


z 


und der Weg y, des ms © Zylinders i 
Ye UAr+ 4) sin ot + (4: — ar cos “N ERBEN: 


4. Diskussion der Verne ee 

Eine besondere Anordnung der Apparatur zeichnete Lissajous-Figuren, die bei einer Phasen- 

verschiebung von 90° zwischen dem angetriebenen und freischwingenden Zylinder in einer 

Gerade ausarteten, die meßtechnisch besonders leicht feststellbar war. Die Konstante von cos at 
in (26) muß in diesem Falle Null werden: 


AA 
Aus (21) bis (24) erhält man die Bestimmungsgleichung für @: 
N YokrÖ j 
(> )-1s = — 026° 4(94,,69)* Tg Pair BP). 2222 2 (AT). 


Vz, ergibt sich aus (19) mit a,,. Die weitere Auswertung erfolgt im En am einfachsten 
graphisch, 

Ist jedoch 9<1, was bei den Danziger Versuchen stets der Fall war, so läßt sich (27) 
rechnerisch weiter behandeln. Man kann näherungsweise setzen: 


Ir, 


te=2 und — P—nal (n=0, 1%. Bnyes: 


Entwickelt man 1er und Tg (Vor öP) in Reihen nach steigenden Potenzen und setzt die 


letzteren Näherungswerte ein, so ergibt sich schließlich für @: 
Or 1+k+0,16K® + --- 


Fe BR : a7 . (28a). 
Darin ist 5 
| 2\2 \ 1 
20,1 (2) a 29a), BEE RE RE ENT 
fi (29a) 9, | (29b) 
und 
„= . eh 14 Gab HOTEL. Due, nr (a0, 


(28a) enthält natürlich ee die Hartungsche Naherıngeläsung: bei der die Eigenträcher 
der Flüssigkeit vernachlässigt wurde, denn mit o—0 geht (28a) über, in Re 


(28b) ergibt für Newtonsche Flüssigkeiten (#—- oo) als kritische Frequenz &,— ®. Die 
Beobachtung einer kritischen Frequenz wurde also als Beweis einer Maxwellschen Flüssigkeit 
angesehen. Bei Beachtung der Flüssigkeitsmasse (0 +0) .ergibt sich dagegen auch für New- 
tonsche Flüssigkeiten eine beobachtbare kritische Frequenz und die Abweichung zur Maxwell- e 
schen Flüssigkeit ergibt nur eine Verschiebung dieser Frequenz. Das Kriterium der Maxwell- i 
schen Flüssigkeit hängt also nicht an dem Auftreten einer kritischen Frequenz, sondern in der 
genauen Erfassung eines (gewöhnlich kleinen) Frequenzunterschiedes. ie Sa 
Für Newtonsche Flüssigkeiten (@ — ) ergibt sich unter Vernachlässigung . der höheren a 
Potenzen von k aus (28a), was-bei den a Versuchen BEL nn ET 


Ar 5 er 
ze el — 0,24 026 ) N 1 
&%,/2isl nach (28b) der Wert, der nach es in Danzig elanhten Ansatz einer 
digkeitsverteilung über den Ringspalt identisch mit fi sein sollte. Bei den a 

ist noch (0,24 0262) neben 1 vernachlässigbar, so daß beiN Ra: sig) 
noch eine Funktion von 0, d und besonders n ist. dee diesem Bauaı 


Eine in Denis entwickelte ee a 
nicht auswertbar. Auch die hier angegebene Verbesserung de: 
kaum etwas ändern. Erst die Verbesserung der an N 
einer gut auswertbaren Amplitudenmethode. RR x: 
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5. Anwendung auf die Danziger Versuche 


Die Genauigkeit der Größe V, ist zur Beurteilung der Flüssigkeit und Bestimmung des 
Schubmoduls entscheidend. Es ist daher notwendig eine Fehlerdiskussion für V, durchzuführen. 
Eine sichere Fehlerdiskussion ist nicht mehr möglich, da man bei der Danziger Näherung die 
Empfindlichkeit der Meßfehler auf die Bestimmung des Schubmoduls nicht sehen konnte. Des- 
halb achtete man auch nicht im notwendigen Maße auf die Kleinhaltung der Meßfehler und 
machte keine genaue Fehlerrechnung. | 

Eine nachträglich gemachte Überlegung, die von Herrn Prof. Fromm und Herrn Dr. 
Hartung stammt, ergab die Möglichkeit eines konstanten Apparatfehlers von -+20%. 
Dazu kommt noch ein Fehler infolge n-Fehlbestimmung, so daß der Gesamtfehler — 25 bzw. 
50% betragen kann. Er dürfte aber kaum über diesen Wert liegen. 

Ist V,—=1 bei fehlerlosen Angaben der Apparät- und Flüssigkeitswerte, so wird @— «, 
d.h. es liegt der Newtonsche Grenzfall vor. Bei den hier möglichen Fehlern kann er auch bei 


; 0,75< V,<s 1,50 
vorhanden sein. In diesem Bereich läßt sich nicht sagen, ob ein Newtonscher, Maxwellscher 
oder gar ein dritter Flüssigkeitstyp vorliegt. 


Ist V,< 0,75, so liefert (28a) einen positiven endlichen Wert für @ und widerspricht nicht 
der Annahme einer Maxwellschen Flüssigkeit. Bei V, >1,50 aber ist ein anderer Flüssigkeits- 


a u 
. > > 


) typus vorauszusetzen. In diesem Sinne sind die noch erhaltenen Danziger Versuchsergebnisse 
ausgewertet worden. Die sich an diesem Abschnitt anschließende Tabelle ergibt eine Zusammen- 
stellung. 

Zusammenstellung der Danziger Versuchsergebnisse. 
m t 10? | @ nach (28b G nach (28 

Flüssigkeit RE I AH ee 

°C Poise g/em? g/em? @ 
Dibuthylphthalat. . 20 0,208 0,020 0,83 ungewiß = 
30 0,132 0,0088 0,92 Ep —_ 
40 0,092 0,0040 0,83 Ri _ 
50 0,066 0,0018 0,73 0,0075 0,04 
“ Dimethylphthalat. . 20 0,175 0,0175 1,06 ungewiß = 
E 30 0,110 0,0067 1,02 AR _ 
40 0,075 0,0026 0,87 fe _ 
50 0,053 0,0010 0,65 0,0029 0,06 
Diäthylphthalat . .| 20 0,130 0,0124, 1,73 Antimax, 
30 0,088 0,0048 1,16 ungewiß — 
. 40 0,065 0,0018 0,81 r _ 
50 0,048 0,0007 0,57 0,0016 0,08 
Glyzerin . . .39% 20 0,036 0,0020 3,85 Antimax, _ 
63% 20 0,137 0,040 6,67 R en 
87% 20 1,45 | 0,80 0,88 ungewiß _ 
3 87% 30 0,77 0,27 0,99 FR _ 
87% 40 0,43 0,095 1,11 5 en 
87% 50 0,26 0,033 1,04 ER _ 

| 87% 60 0,17 0,0115 0,86 AN er 

DES 2; 20 6,52 7,50 0,27 10 0,2! 

% 30 2,95 2,50 0,45 4,5 0,13! 

40 1,42 0,85 0,67 3,6 0,05 

50 0,75 70,30 0,83 ungewiß == 

60 0,42 0,10 0,89 ER _ 

a RE 40 14 22 0,17 36 0,2! 
50 4,7 6,5 0,47 12 0,1! 

Au | 60 1,82 2,0 0,93 ungewiß | _ 


F Man ersieht aus ihr, daß V,in der Mehrzahl zwischen 0,75 und 1,50 liegt und man hier 
keine Angaben über den Schubmodul machen kann, ja selbst der Flüssigkeitstypus ungewiß ist. 

_ Fünf Flüssigkeitsarten (acht Angaben) können Maxwellscher Natur sein. In diesem Falle 
_ dürften die hier gemachten Angaben über den Schubmodul richtig sein. Bei Diäthylphthalat und 
en 39% und 63% muß dagegen bei 20°C ein dritter noch unbekannter Flüssigkeitstypus 
- Interessant ist, daß Diäthylphthalat bei 20°C dem Maxwellschen Stoffgesetz wider- 
Dieser Widerspruch nimmt mit steigender Temperatur immer mehr ab und bei 50°C 
ie Flüssigkeit von Maxwelltypus zu sein. Es könnten daher in einer Flüssigkeit 


88 Pieruschka, Die Grundlagen einer Meßmethode des Schubmoduls ne ee 3 a u 


Eigenschaften vom Maxwellschen und Antimaxwellschen Charakter vereint sein, die beide in 
verschiedener Weise von der Temperatur abhängen. So könnte einmal der eine und das zweitemal 
der andere Typus überwiegen. Ähnliches scheint bei Glyzerin vorzuliegen. 


III. Abschätzung der durch die Linearisierung entstehenden Fehler 
6. Verbesserung der linearisierten Lösung 


Unter der Voraussetzung, daß die Amplitude der ersten Oberschwingung klein gegen die 
der Grundschwingung ist, kann die Lösung des linearisierten Systems als eine erste bereits recht 
gute Näherung für 9v/dx (25) angesehen werden. Berechnet man diese für die Stelle des frei- 
schwingenden Zylinders und führt sie in das Grundsystem (3) und (4) ein, so erhält man ein 
System gewöhnlicher linearer Differentialgleichungen für o und tr. Die Randbedingung (6) 
liefert eine einfache Gleichung zur Berechnung der Beschleunigung 5b des freischwingenden 
Zylinders. Die Gl. (5) ist hierbei dadurch berücksichtigt, daß sie zu dem Ergebnis (25) bei- 
getragen hat, welches jetzt in (3) und (4) eingesetzt wird. 

Für die Verhältnisse am freischwingenden Zylinder (e—0) läßt sich (25) unter Verwendung 
von (16) und (17) in der Form schreiben: 


5 
rn Fasinat, ei), 


wenn 

e EN FHNlET. 

AT RR | 

u 3 *__ 4* * * 

E RR at SI ARE. 70 On Bine Or (32) 
4 (43 +4) — 9 (45 — A}) 

2 & 
Ri; ae | | 

z a h, gesetzt wird. Nennt man die Beschleunigung am freischwingenden Zylinder b, so folgt aus (6) 


t=ub. Setzt man letzteren Wert und (31) in (3) und (4) ein, so erhält man das simultane 
lineare Differentialgleichungssystem 


; 1 

+ =, in a ee = Ge SIND A el 
N KEREHNR 6 \ | 

+7 '0=—yuabsinat, 2 ME En EEE ach 


. 


TIER R 
vr an 


Dieses läßt sich nach dem D’Alembert-Ampereschen Verfahren #) Jösen. N 

Wendet man diese Integrationsmethode auf den hier allein interessierenden Burn 3 
zustand, bei dem die Ausdrücke mit den beiden Integrationskonstanten Null werden, an, so 
Ei man für b und oe: re 


be 7 (cos A: ertl? [ eıT sin Adt, — sin A eh? [ eilt cos A dt) . 2% (35),. 


Kr | 1 a 
= — el sin Aettl? | eu? sin Adt, 4 q cos A ‚tlR font dı,) (36). RE 


Hier ist 2— cos a. | ER 

7. Die Beschleunigung des freien Zylinders 

Um eine für die numerische Rechnung für die Beschleunigung b brauchbare Fo 
War löst man cos und sin A nach steigenden Potenzen von h auf jetzt, 


durch 4 Fi cos &t;j. Dann werden die ganzzahligen Potenzen von cos atı nach 


Vielfachen von at, entwickelt. Auf diese Weise erhält man: 


", 
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r 


mit 


n—=N, % 2, 3 
27-4 7\ (a1/2 S)?iH ; =0undl, 
E* —p, Hu See Be 1 
2n-+j 7, Pr. ZN ) i—N -AiHjt j med für ne12 
n=1. 


Setzt man die zuletzt angegebenen Formeln, in die Integrale von o) ein und trifft die 
Abkürzungen | 


ae Re 
9; — IrG9% ir A (37a), und Ti Tears 


so gewinnt nach einer leichten Integration (35) die Form: 


» 


* 


. (37b), 


a; 


/ male SEE una) Team+ı) cos 2 nat, sin (2m +1) at, 
4 \n=0 m=0 


- 223 ee Tom: cos (On 1) ot, sin 2mat, 
m = 


Ss 
IR 

o 
o 


N RER RE ESTER Un eh 


+2 I Bi na Oam+ı)cos2 nat, cos (2m -+ 1)atı 
n=0 m= 


N 3 — Eim Eeen+n O;m - cos 2m at, * cos En+1)ahl. 
r n=0 m=0 
Nach Vertauschen der Reihenfolgen der Summation und Zusammenfassen erhält man, 
wenn man noch die Produkte der trigonometrischen Funktionen durch die Summen ersetzt: 


a Fol 
& 


Saul 53 Ein Eim+1 (Tamtı a Tan) sin [2 (m +n)+ lat, 


n=0 m=0 


en a en 


sn 3 E2n Eim+ı (Tonsı + T3„)sin [2 (m Er 1]atı 


n=0 m=0 


+2 2 Ben Bim+ı (Iam+ı — Opn) cos [2 (m + n) +1] otı 


2 | 

A a z 3 „Bin Ei m+1ı (O2 m+ı — an) cos [2 (m — n) + 1] hl. 

.’ { n=0 m= 

2  Esist zweckmäßigin den Ausdrücken Deifın [2 (m—n)-+1] at, und cos [2 (m —n)-+-1]at, 


> ‚das Vorzeichen von n zu wechseln und dann über » von — oo bis 0 zu summieren. Begantet man 
; noch, daß 

"1 Eu i “ ur Diana — Ion ‚und RATEN 

ist, setzt weiter formal 


BR ae N a a Te Dan Kur 4028 
fest und nimmt N 
ei ne) 41/2 6)?°+5 .n=0,1,23,... 
Band Eu u PRREET ud ü a 


E3; ni so erhält man 
le als 2 Bantı Bun (Tamıı — Tan)sin 2m tn) + Nat, 
x Mm=UN=— oo 


RR Be E Bantı Ban (Oumsı — Bun) Bm tm + al. 
m= —w 
man noch N Be - Mm bzw. N BR —_ N — m, so erhält man nach ein paar ein- 


Eau sin @ k + 1). at + 2, Bascı [KoR} (2 k + 1) sul . B (40) 
47 e ? Em x ae et. 
N NE R 
dr Banır eins Er Tenn)—E, m+1Bam+2x+2 (Tamt+ı +Tam+2x+2)} (41) 


2% 
5 SE RN “ 


ar: 


EN a RAR 
ee“ Bade 7) 0 ‚ y Berler: ne 


i x a e: PEN , ra j r3 5 Rn) a \ 
; 2 2 
90 Pieruschka, Die Grundlagen einer Meßmethode des Schubmoduls „7 angew, Math. Mech. 


Bd. 31 Nr. 3 März 1951 
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Damit ie das Problem gelöst. Mittels (37), (38), (39), (41) und (42) kann man nach Eon b 
berechnen, durch Integration von diesem v, und durch ng Integration %9- 
8. Die Normalspannungen am freien Zylinder 


a) Berechnung 
co berechnet sich am einfachsten, wenn man mit dem Ansatz 


= GI (Drrsin 2 kat, + Orr cos 2katı) re ee be 
k=0 
und (40) in (33) und (34) geht. Man erhält ein lineares Gleichungssystem, aus dem man die Dy;. 


und Cs; berechnen kann, wenn die Ag;ıı und Ba;+ı bekannt sind, Von Interesse sind nur 
die ersten drei Werte. Man erhält für diese: 


Kit 57 12 Eee 7} 
29 4 
(,= +2 (Brom, —B)— (TA) ae 
29 
Dt ml) 9a AL (BIER 


Die Normalspannungen lassen sich auch unmittelbar aus der Lösung der linearisierten Form 
berechnen. Im allgemeinen werden diese Werte nicht mit (43) übereinstimmen. 


b) Anwendung auf die Danziger Versuche 


Theoretische Erwägungen und numerische Rechnungen zeigen, daß bei nicht zu großen 
Ö-Werten sich das nach (32) berechnete a, nicht nennenswert von der Amplitude «a des antrei- 
benden Zylinders unterscheidet. Daher soll in diesem Abschnitt a, gleich @ gesetzt werden. 
Weiter zeigt eine einfache Überlegung, daß diese Näherung bis auf die Phasenverschiebung 
identisch mit der Annahme eines konstanten Geschwindigkeitsgefälles im Ringspalt ist, wenn 
die Amplitude des freien Zylinders etwa unter 10% von der des angetriebenen Hleibt, was bei 
den Danziger Versuchen erfüllt war. . 

Mit dem Wert der Danziger Apparatur a/2 ö = 3,67 erhält man für die Werte der #, 
nach (39): 


B, u a a 


0,2847 0,0935 0,2592 0,2348 0,0673 0,3081 —0,3525 | —0,2682 
L; E, En Eu 
| 0,1580 0,0785 | —0,0334 | —0,0126 


12 E 13 E 14 E 15 
0,0013 —0,0004 | —0,0001 


' Nun sind nach (41) und (42) A, 
As B,, B; berechnet und in Bild2 
als Funktion von ® dargestellt. Da- 
beiistauch die Rechnung für größere 
- 9-Werte durchgeführt, als es bei den 
getroffenen Näherungen zulässig ist. 
Das Charakteristische bleibt aber auch Fe 
da erhalten. 
Man erkennt, daß die Amplitude 
der Grundschwingung praktisch durch 
A, wiedergegeben wird. B, hat nur a 
auf die Phasenverschiebung® einen IE 
Einfluß, die hier von re beson- 
deren Interesse ist. 


Bild 2 


Berechnung von. 
und ul nach 
2 s inte 
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r FR dasselbe Bild ist das prozentuale Verhältnis der Weg-Amplitude der ersten Oberschwin- 
i gung zur Grundschwingung aufgenommen. Es ist mit A, bezeichnet: 
Y4:+ B: 
mehr: Fr u BEE Le (47). 
Y4+B: 


Man erkennt aus der Darstellung, daß A, nur wenige Prozente beträgt. Die Oberschwin- 
gungen treten neben der Grundschwingung nicht sehr in Erscheinung. 


c) Kritik der Näherung o/@<1 

Zur Beurteilung des Zulässigkeitsbereiches der hier getroffenen Näherung ist der Wert o/G@ 

von Bedeutung. Bei der Linearisierung des Differentialgleichungssystems im Abschn. Il ist o/@ 
neben 1 vernachlässigt worden. Über die dadurch entstehende Ungenauigkeit gibt Bild 4 einen 
Anhalt. Die Amplitude der mit der Kreisfrequenz 2 wechselnden Normalspannung ist gleich 
dem konstanten Anteil dieser Spannung. Daher kann man bis 
zu einem gewissen Betrag von o/@(1-+0/@) durch die Kon- 
stante (1-+ ou(@)/@) ersetzen. Dabei ist aber zu beachten, daß 
o,(2) noch eine Funktion von z ist. Man kann annehmen, daß 
für |o,/@| < 0,2 die Rechnung befriedigendes Vertrauen haben 
kann. Nach Bild3 ist 90,1 für |o/@| = 0,2. Demnach gilt 

unsere Rechnung für 


<0l. Bilde 


Diese Bedingung war bei den Danziger Versuchen fast stets erfüllt oder wurde nur wenig über- 
schritten, so daß die Linearisierung hier zulässig ist. 

Man kann noch höhere d-Werte zulassen, wenn man in (3) (1-+0/@) durch (1-# ogmitier|F) 
ersetzt. Dabei bedeutet oy„.r.; einen Mittelwert des konstanten Anteil von a —0o,(#) — über 
die Breite des Ringspalt genommen. Überschlagmäßig kann man o9„..1/@ durch o,/@ ersetzen. 
Dann können die Ergebnisse der Rechnung des linearisierten Systems (Abschn. II) übernommen 
werden, wenn man für @ einen konstanten Abminderungsfaktor einführt. Setzt man also in (3) 
statt @G @G, mit 


Ga lA88), | m—N']ır We (45h) und J=1+o/l@ .- - (48c), 
so wird (3) zu 
aa „2,80 
Gasen. az 


Aus der letzten Form der Gl. (3) erkennt man, daß man eine gewisse Berücksichtigung 
von 1--0/@ erhält, wenn man in den Ergebnissen des Abschn. II @ durch @, und n durch n, 
ersetzt. Man wird aber auch hier eine gewisse Einschränkung für o,/@, d.h. f vorschreiben 
\ müssen, die ungefähr bei 
- \oo/@] S 0,4 
angesetzt werden kann. Nach Bild 3 kann man setzen 
ee DR ee eek in 3 (ASc) 


mit der Bedingung 9 < 0,20. Damit wird (28a) zu 


(28d) Er A mit (30b) yo (V, nach (30a)). 


| d-n)f f 
ö 28d) ist durch Iteration zu lösen. 
® Aus der Form von (30b) ersieht man, daß die negativen Werte für @, die sich aus den 
E ‚Danziger Versuchen bei Anwendung der Formel (28a) ergeben, nicht durch die verbesserte 


Näherung nach (28d) positiv werden können. Da f <1 ist, wird V nach (30b) bei von 1 merkbar 
Fin a seniedlichen Werten von f noch größer als V,. 

» 9 dung der Amplitudenmethode 
$, _ Im Abschn. 5 wurde erwähnt, daß eine Auswertung der Amplitudenmethode tele der 
ıc echnung für das linearisierte System an der für die Streuung der Meßwerte zu große Empfind- 
ichkei dieser Bahhde scheitert. Von entseheitlendem ‚Nachteil ist hierbei, daß die Lineari- 


setten. Rechnung (Abschn. III) an man auch größere Deere für 
tung zu lassen und adarch eine Methode nicht nur zur Abschätzung des Schubmoduls, 


- 

nr 

> 
22 
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sondern insbesondere zur Kontrolle für die Zulässigkeit des von Fromm erweiterten Maxwell- 
schen Stoffgesetzes gewinnen. 

Zur Auswertung der über die Kreisfrequenz aufgetragenen Amplitude einer Meßreihe geht 
man vom Abschn. 9b) aus. Aus den Bildern 2 und 3 ersieht man, daß die Amplitude im wesent- 
lichen durch A, bestimmt wird. Unter dieser Voraussetzung erhält man nach (40) für die an 
tude A, des freien Zylinders: 


a BEL (49a) ‚ oder GA Et a5 N Fein (49b). 


Man kann also nach (49b) @ A, über die Kreisfrequenz auftragen, da u, « und A, bekannt 
oder meßbar sind. Durch Vergleich dieser graphischen Darstellung mit der berechneten Kurve 
A,(®) läßt sich @ leicht ermitteln. Am einfachsten wird diese Berechnung von @, wenn man 
vom Extremwert der gemessenen A,-Kurve ausgeht. Bei voller Zulässigkeit der Rechnung ist 
Amar = 0,295 und liegt bei « = 
nung von @ und kann dadurch die Rechnung gleichzeitig kontrollieren. Nennt man diese aus- 
gezeichnete Kreisfrequenz &,, so lauten die Formeln zur Berechnung von @: 


A, ; a 5 
0,59 bzw. G=6n%, u oe ee ee Io re € 


Nun kann man A, als Funktion von ® darstellen. Eine Übereinstimmung mit der berech- 
neten A,-Kurve würde die Berechtigung des Verfahrens beweisen. Sollte dieser Vergleich nicht 
befriedigend sein, so kann man einen Schritt weiter gehen und untersuchen, wieweit man die 


PAST 


— 0,17. Damit hat man gleich zwei Gleichungen zur Berech- 


G= 


02 


01 


0 05 10 15 20 25 


7) DR 10%. 15 20 25 y i 
Bild 5 ; Bild 6 e 


Versuchskurve durch“ein Abnehmen des Verhältnisses = mit steigendem 9 (»>0-% 2) 


erklären kann. Eine Übersicht hierzu bieten die Bilder 5 und 6. Hier sind A, und u: 

über 9 für a,/a= 1, 1/2 und 1/10 aufgetragen. Die Kurven gelten für die Danziger Versuchs- - 
anordnung. Man erkennt an Hand des Bildes 5, daß das Maximum sich mit abnehmenden a,fa 9% 
in Richtung größerer d-Werte verschiebt und dabei nur wenig abnimmt. Aus Bild 6ersieht man, 4 
daß mit abnehmenden a,/a auch, —o,/@ stark abnimmt und sich dabei der ren 
der Näherungsrechnung erweitert.‘ | 
Eingegangen: 25.4. 1950. N | . BRRE: f I. 
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Ein. hinreichendes Kohvergenzkritörkum und Bei der Iteration in Gesamtschritten?) rechnet. 


eine Fehlerabschätzung für die Iterzitlon in man die Folge der end a für die wi 
Einzelschritten bei linearen Gleichungen.) folgt: wi hal 


Wir betrachten das System linearer ae n BARS 1] 

+) _ —_t 
a1%y 4 Apt2t+ "+ ann + a, = 0 “ ai | Q 
a21%, 7 A22Xa tr en E21) et mn 
EEE N ae a IR SE SR a >= 2 ER 
anı&ı + Angta+ "+ amnzn + an= 0 


RE ir BEER 3 Beim Einzelschrittverfahren = 
mit den Un nnten ©, bis ©, und der Deter ninante 
det a;z-+0. Für das Iterationsverfahren sollen neue, ande die bereits a 
für die Praxis nützliche hinreichende Konvrergenz- 
kriterien gegeben werden. 


’) Aus dem Institut für Praktische Mathematik 
nischen Hochschule Darmstadt, 
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nach unten“ iteriert, ergibt sich die Vorschrift 


er el ve (v-+1) Ye 
u = |Zaum + 3 ou + 
ü l=1 (k=i+1) 
Geht man von der n-ten Gleichung aus, iteriert also 
„von unten nach oben‘, so hat man in (B) nur auf 
der rechten Seite Bir + mit a) zu vertauschen, 


Für die Verfahren (A) und (B) stammt von Collatz?) 
die gemeinsame hinreichende Konvergenzbedingung 


(B). 


=u<1l...0) 


n 
Max u; = Max N 
v % k=1 
j ki 
mit der Fehlerabschätzung 
w41) 21 [ii 
Max ®, a" %; = 1 way 


Der im allgemeinen besseren Konvergenz des Einzel- 
schrittverfahrens (B) werden die folgenden hinreichen- 
den Konvergenzkriterien mehr gerecht: 


air 
ii 


Max 125” = an (3). 
® 


' Kriterium I. Das Einzelschrittverfahren (B) kon- 


vergiert, wenn 
Mxu;=a<l...... ee (&) 
i 


bleibt, wobei die &; durch die Rekursionsformel 


Ba 
He 


5 n 
Teil Zlaul &% +, 2 laml De) 
erklärt sind. Dabei gilt die Fehlerabschätzung (3), 
wenn dort u durch a ersetzt wird: 


& 
1—& 


Max | +" _ «| = Max [2 ar |» (6). 
ı ® 


Die &; werden also nach der Iterationsvorschrift (B) 
unter Fortlassen der ai, Nehmen der Absolutbeträge 
und mit Ausgangswerten 1 berechnet. 


Kriterium II. Das Einzelschrittverfahren (B) 
konvergiert, wenn | 
Maxß;=Pß <1 
® \ 


bleibt, wobei die f; durch die Rekursionsformel 


(Zu) re 3 |air | . (8) 


wi k=i+1 
‚erklärt sind. Dabei gilt die Fehlerabschätzung (9) 


| 1 
a lei: | 


ee ß 
letz 
Die ß; entstehen also durch Vergröberung der a;. 

indem man während der Rekursion mit dem Maxi- 
mum der bereits gefundenen ß, weiterrechnet; da- 

durch lassen sich vorbereitend die Koeffizientenbeträge 
‚links und rechts der Hauptdiagonale summieren und 


Max jet? _ «| (9%. 
® e 


Be somit Multiplikationen einsparen. 


. Eine weitere Vergröberung würde darin bestehen, 


0 in (8) rechts Max ß. = 1 zu nehmen; dadurch kommt 


# 


man auf die bisherige Bedingung (2) zurück. Die 
Verbesserung 


durch die neuen Kriterien zeigt sich 
besonders deutlich am Falle eines gestaffelten Systems, 
_ bei dem rechts der Hauptdiagonale alle a0 sind. 


tn. 


ben nämlich die Kriterien I und II die Werte 
05 . sie lassen erkennen, daß in diesem 
Iterationsvorschrift die exakte Lösung in 
‚Schritt liefert. Das alte Kriterium (2) zeigt 
die exakte Lösung erst dann an, wenn alle 


bschätzung für das Iterationsverfahren en vr i 
_ das Kriterium (2). 


TUR 
rt 
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ar) = + ye SET (11) 
folgt aus (B) und (1) 
BE x 
eu: a Dar "+. 3 0% N . (12) 
ii Ik-1 k=i-+1 
und weiter (13). 
a en un Sn aaleleol 
In! | Es a; u jGik, /%% Kr P> |aik| h wi 
“le k=i-+1 
Setzen wir 


Max |y)|=& und 5 Lern (14), 
i 


so wird 
1 i—1 n { 
d; Sl Zlenla+ P} ja . (15) 
lei | 15T k=i+1 
Ein Vergleich mit der Folge (5) lehrt 


6;<oa; und re 6< le. I LE) 
Da nun aus (14) folgt 
Ds ve+® = er RAR REN (17), 
besteht nach (16) die Ungleichung 
2 urn) <ae=a nn WM JE 


also 


Max wen] <o”*1Max |. 2 (10E 
ı (7 


Hiernach ist für a<1 Konvergenz gewährleistet. 
Zur Fehlerabschätzung bilden wir 


Max |.” #2 a’ a 8 a. Max |yy”) en rn | 
> Max ve] — Max yer+D] | 

und finden mit (18) 
Max a +n _ 2m) En Max (+2) | 


— Max + an 
ENT 


(20) 


(21) 


Es gilt also 
Max at’ Watt “| = Max |yi’* S 
ı ı 
.. (22) 
% w+1) _ „o)) 
< ar Da et | 


oder (6), was zu beweisen war. 


Zahlenbeispiel. Wir schreiben zur besseren 


Übersicht nur das Schema der Koeffizienten an und 
nehmen als Ausgangswerte aD für die Iteration nach 
‘dem Einzelschrittverfahren (B) die a;: 


ar. re 


Tatsächlich sind die 2, 


. 
. 
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Anwendung des Kriteriums I ergibt 


= u 0,2+]-0,1140]= 03, 

AR Br [0,8-0,3 + 0,14 | 0,11] = 0,44, 

na - [I 0,9]-0,3 + 0,1-0,44+ 0,1] = 0,414, 
= 1): 0,3 + | 0,3]-0;44 + 0,4 0,414] 


—= 0,4776, 


Die Feststellung Maxa;=a » 0,48 <1 erweist 
das Verfahren (B) als konvergent. 

Das Kriterium II liefert a,=f, wd ,=ß,, 
weiter 


Ba= SE [I— 0,9] + 0,1) 0,44 + 0,1] = 0,54, 


Wegen Max Pi=Bß= 0,702 <1 zeigt auch das 


gröbere RAterinen II noch die Konvergenz des Ver- 
fahrens (B) beim vorliegenden nein an. 
Für die Fehlerabschätzung ist (3) überhaupt nicht 
anwendbar, (6) liefert für 29) 


2) & (1) 


Ds SE A Ka) et 
Er ®; x ri ®; 
ze 0,17<£0,16. 


477 
;”, wie ihr Vergleich mit den 
exakten x; zeigt, genauer. 


Darmstadt. H. Sassenfeld. 


Zur elementaren Theorie der Zapfenreibung. 


In vielen Lehrbüchern über technische Mechanik 
und technischen Handbüchern wiederholt sich bei der 
Behandlung der Zapfenreibung in allen Auflagen ein 
anscheinend bisher merkwürdigerweise nicht bemerk- 
ter Fehler. 

Der Zapfen vom Radius r wurde durch die gegebene 
Kraft P gegen das Lager gepreßt. Das Moment M, das 
nötig ıst um den Zapfen zu drehen, bzw. dasim Grenz- 
fall des Haftens gerade noch Gleichgewicht hält, wird 
angesetzt: M = r&4AR, wobei die Reibungskräfte 
AK gleich dem Produkt aus der Reibungszahl u und 
den Normaldruckkräften AN sind, deren Verteilung 
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jedoch im allgemeinen unbestimmt ist. Um wenig- 


stens eine untere Grenze für M festzulegen, wird nun 


in der in Frage stehenden Ableitung das Gleichgewicht 
des ruhenden Zapfens ohne Reibung herangezogen, 
für welches gilt: 


.weil die algebraische Summe der Seiten or Kräfte- 


polygons der AN größer ist als seine Schlußlinie ?. 
Damit kann geschrieben werden: 


Mi aut Pin su (2), 
wobei für die Zahl a gilt: 
NH ur Er (3). 


Die Anwendung der Gleichgewichtsungleichung (1) 
auf den Fall mit Reibung ist unzulässig. Für diesen 
Fall gilt vielmehr analog: 24K>P, wo AK die 
Resultanten aus Normaldrucken und Reibungskräften 
bedeuten (zur Gleichgewichtsbetrachtung denke man 
sich P mit M zu einer gleich großen gleichgerichteten 
Kraft zusammengefaßt, deren Wirkungslinie nach 
rechts verschoben ist). Mit Einführung des ar; 
winkels e durch die bekannte Definition: tge = 
wird AR=AKsino und M ergibt sich zu u 
rPsino. Wenn daher der Ausdruck für M in der 
Form (2) beibehalten werden soll, ist Ungl. ( e) richtig 
zu ersetzen durch: 

& >cosp. 


Der Grenzfall x = cos o (der nach (3) unmöglich 
wäre) tritt ein, wenn in einem Lager mit Spiel die Be-, 


rührung nur längs einer ‚Bizpugenden des Zapfens er- | 


folgt.t) 


Tucumän. E. Mönch. 


3) Der gleiche Gedanke ist schon in dem Buche H.Blasius 
Mechanik I, Hamburg 1948, ‚geäußert worden, _ W; 
j 


BUCHBESPRECHUNGEN 


Dr. Ing. R. Zurmühl, Matrizen. Eine Dar- 
stellung für Ingenieure. ’xv -+ 427 Seiten mit 25 Abb. 
Berlin-Göttingen-Heidelberg 1950. Springer-Verlag. 
Preis geb. 25,50 DM. 

Matrizenrechnung ist vielen Ingenieuren heute noch 
unbekannt, oft sogar auch dem Namen nach, und 
doch verdient sie ebenso beachtet zu werden wie an- 
dere große Gebiete der Mathematik, z. B. die Vektor- 
rechnung. Die Matrizen bilden ein übersichtliches und 
weitreichendes Hilfsmittel für alle Untersuchungen, 
bei denen mathematisch gesprochen lineare Systeme 
zu behandeln sind, z. B. Systemelinearer Gleichungen, 
linearer Differentialgleichungen, Systeme von Strom- 
und Spannungsgrößen in elektrischen Netzen usw. 
Auch z.B. bei kleinen Schwingungen sind die Matrizen 
das dem Problem angepaßte Mittel, und es wird wohl 
gerade bei den mechanischen "und elektrischen 


Schwingungen die Überlegenheit der Matrizen dazu 


auch im Hochschulunterricht den ihnen grliährenden FR 

Platz einnehmen. 
Es gibt über Matrizen bis. din nur ler o lehr- 
ER 


rung für ne 
führlich und weitzeichend 
zierteren Anwendu 2 


Durch diese A 


führen, daß sie in der Lehrbuchliteratür und. dann 
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Mathematiker, sondern will den Ingenieur mit dem 
neuen Gebiet vertraut machen, ohne ihn mit langen 
strengen Beweisen zu ermüden. 

Das I. Kapitel bringt Grundlagen der Matrizen- 
algebra, das II. Kapitel Determinantensätze, Rang 
einer Matrix und lineare Gleichungssysteme, das 
III. Kapitel Formen und linesre Transformationen, 
das IV. Kapitel das Eigenwertproblem, erst für reelle 
symmetrische Matrizen, dann für beliebige komplexe, 
mit „speziellen“ und ‚allgemeinen‘ Eigenwertauf- 
gaben, und Paare quadratischer Formen, das V. Ka- 
pitel die Struktur einer Matrix, Elementarteiler- 
theorie, Zurückführung auf verschiedene Normal- 
formen, Hauptvektoren, Matrizenfunktionen. Das 
VI. Kapitel bringt von den zahlreichen bisher auf- 
gestellten numerischen Verfahren die. wichtigsten 
Methoden zur Lösung linearer Gleıchungssysteme und 
Bestimmung der charakteristischen Zahlen von Ma- 
trizen (bei Schwingungsaufgaben ist das die Bestim- 
mung der Eigenfrequenzen). Das VII. Kapitel ist den 
Anwendungen gewidmet, davon der erste Paragraph 
der Elektrotechnik. Der nächste sehr wichtige Para- 
graph gibt eine. geschlossene Theorie der kleinen 
Schwingungen. Es folgen Systeme von Differential- 
gleichungen mit konstanten Koeffizienten, wieder mit 
} zahlreichen Zahlenbeispielen, Differentialmatrizen, 

nichtlineare Transformationen, Tensoren und Ma- 
trizen in der Ausgleichsrechnung, Vielleicht ist manch- 
mal der Verfasser in seiner Liebe zu den Matrizen 
etwas weit gegangen, z.B. bei der Darstellung der 
Tensoren, bei denen der Referent lieber den in der 
Mathematik entwickelten Tensorkalkül verwenden 
würde, aber der Verfasser hat sich ein großes Ver- 
dienst erworben, indem er in dem ganzen sich auf die 
Anwendungen beziehenden Kapitel häufig Dinge, die 
schon bekannt waren, durch Umschreibung in die 
Matrizensprache in neuer einfacherer Weise darstellt. 

Man möchte dem sorgfältigen und gründlichen 
Buche wünschen, daß es von recht vielen Ingenieuren 
benutzt und möglichst auch durchgearbeitet wird. 


Hannover. Collatz. 


Dr.-Ing., Dr.-Ing. e. h. R. Grammel (Prof. a. d. 
Techn. Hochschule Stuttgart), Der Kreisel, 
seine TheorieundAnwendung. Zweite 
neubearbeitete Aufl. Erster Band: Die Theorie 
desKreisels. XI-+ 281S. m. 137 Abb. Berlin- 
Göttingen-Heidelberg 1950. Preis brosch. 30,— DM, 
geb. 33,— DM. 

Das Grammelsche Buch über den Kreisel war seit 

vielen Jahren vergriffen und sehr gesucht. Allgemein 

) wird daher sein Wiedererscheinen begrüßt werden, 
- umso mehr, als es durch eine vollständige Umarbeitung 

ganz wesentlich . gegenüber der alten Auflage ge- 

wonnen hat. Die beiden Abschnitte der früheren Auf- 


lage sind jetzt zu je einem Bande erweitert. Davon 


liegt der erste, der die Theorie bringt, vor. Er beginnt 

mit einer ausführlichen Darstellung der für die Kreisel- 

| theorie nötigen Grundlagen der Vektorrechnung und 
r der Mechanik und mit einer Ableitung der wichtigsten 
‘ Eigenschaften des Trägheitstensors, bringt dann im 
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frei, geführt, unter Zwang und Stoß und schwer) und 
anschließend im dritten den kräftefreien und den 
schweren unsymmetrischen Kreisel. Hier wird auch 
die Staudesche Drehung sowie dıe Nachbarbewegungen 
der stabilen dieser Drehungen behandelt. Der vierte 
Abschnitt bespricht interessante Fälle des Kreisels 
im weiteren Sinne wie Spielkreisel, Gyrostat, rollendes 
Rad und rollende Kugel, behandelt ferner die Theorie 
allgemeiner gyrostatischer Systeme. In einem Anhang 
endlich wird die Kreiselbewegung durch Theta- 
funktionen dargestellt. So ist ein didaktisch außer- 
ordentlich geschickter Aufbau der Kreiseltheorie ent- 
standen, bei dem immer der Hauptwert darauf gelegt 
ist, einen Einblick in die mechanischen Vorgänge zu 
gewinnen und bei dem die mathematischen Ent- 
wicklungen überall diesem Ziel untergeordnet sind. 
Auch im einzelnen kommt überall die reiche Lehr- 
erfahrung des Verfassers zur Geltung. Dadurch, daß 
die Ergebnisse der Entwicklung jeweils in kurzen 
prägnanten Sätzen zusammengefaßt sind, wird die 
Durcharbeitung des schönen Buches wesentlich er- 
leichtert. Erfreulich ist ferner, das noch in diesem 
Jahre auch der zweite Band, der die Anwendungen. 
bringt, erscheinen wird. 


Dresden. Willers. 


H. Beer, W. Flügge, K. Karas, K. Marguerre, Th. 
Pöschl, E. Reißner, A. Schleusner, S. Timoshenko, 
E. Chwalla, 0. Dediö, E. Friedrich, 0. K. Fröhlich, 
P. Funk u. E. Berger, 6. Heinrich, K. Kriso, A. Leon 
u. F. Vitovee, W. Mudrak, F. Müller-Magyari, H. Par- 
kus, A. Pucher, F. Reinitzhuber, A. Slibar u. F. Vito- 
vec, F. Söchting, E. Tschech, E. Tungl, H. Winter. 
Beiträge zur angewandten Mechanik. 413 8. 
mit zahlreichen Abb. Wien 1950. Franz Deutike Ver- 
lag. Preis: 30,— DM. 

Aus Anlaß des 65. Geburtstages des Ordinarius für 
Allgemeine Mechanik, Hydro- und Aeromechanik an 
der Techn. Hochschule Graz Prof. Dr. Federhofer 
und des 60. Geburtstages des Ordinarius für Tech- 
nische Mechanik an der Techn. Hochschule Wien Prof. 
Dr. Girkmann haben ihnen Freunde und Schüler 
eine Festschrift gewidmet. Aus zeitbedingten Gründen 
wurden die Arbeiten auf photomechanischem Wege 
vervielfältigt; auch so ist ein außerordentlich sauberes 
und übersichtliches Satzbild entstanden. An den An- 
fang des Bandes sind kurze Lebensläufe und Schriften- 
verzeichnisse der beiden Jubilare gesetzt. Es folgen 
26 Abhandlungen, die sich in der Hauptsache mit 
Problemen aus den Arbeitsgebieten der beiden Herren 
beschäftigen. Man findet sechs Arbeiten aus dem 
Gebiet der Baukonstruktionslehre, sechs, die sich mit 
der Theorie der Platten und Schalen beschäftigen, 
sieben behandeln andere Probleme der Elastizitäts- 
theorie, drei Eigenwertprobleme und die restlichen 
vier verschiedene Fragen der Kreiseltheorie, der Boden- 
mechanik usw. Auf die einzelnen Beiträge einzugehen, 
fehlt hier der Raum. Doch sei betont, daß der vor- 
liegende Band eine Reihe sehr wertvoller Abhand- 
lungen bringt. 


| zweiten Abschnitt den symmetrischen Kreisel (kräfte- Dresden. Willers, 
1 Ü Die besprochenen und angezeigten Bücher sind durch den Buchhandel zu beziehen. 
Ei 
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2 > ; ‚Bei der Schriftleitung sind folgende Bücher eingegangen (ausführliche Besprechung bleibt vorbehalten) 


Me  H. Heyde (Prof. a. d. Humboldt-Univ. Berlin), Me- 
 ehanikfür Ingenieure. Bd. 1: Statik/Dyna- 
milk (Teubners Fachbücher für Maschinenbau und 


» 


100, 


'eubner, Verlagsgesellschaft. Preis geb, 
„ = Fi ri ei L " y) 


hnik). X -+ 322 $. mit 317 Abb. Leipzig. 


K.Federhofer (o. Prof. a. d. Techn. Hochschule 
Graz, Prüfungs- und Übungsaufgaben 
aus der Mechanik des Punktes und 
des starren Körpers. Teill: Statik. 
V -+ 130 S. mit 243 Textabb. Wien 1950. Springer- 
Verlag. Preis 9,60 DM. 
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Dr. F. Neiss (Prof. a. d. Humboldt-Univ. Berlin), 
Analytische Geometrie. VIII + 167 S. 
mit 64 Abb. Berlin - Göttingen-Heidelberg 1950. 
Springer-Verlag. Preis brosch, 9,60 DM. 


6. Pickert, Einführung in die höhere 
Algebra (Studia Mathematica/MathematischeLehr- 
bücher Bd. VII). 298 S. Göttingen 1951. Vanden- 
hock & Ruprecht. Preis brosch. 12,80 DM, geb. 
14,80 DM. 


Dr.-Ing., Dr.-Ing. e.h. R. Grammel (Prof. a. d. TH 
Stuttgart), Der Kreisel, seine Theorie 
und seine Anwendungen. Zweite neubear- 
beitete Auflage. Zweiter Band. Die Anwen- 
dungen des Kreisels. VI + 268 8. mit 
133 Abb. Berlin-Göttingen-Heidelberg 1950. PREIS 
Verlag. Preis geh. 30,— DM, geb. 33,— DM 


Dr.-Ing. habil. Kurt Nesselmann, Die Grund- 
lagenderangewandten Thermodyna- 
mik. XT + 320 8. mit 311 Abb. und 5 Diagrammen. 
Berlin-Göttingen-Heidelberg 1950. Springer-Verlag. 
Preis geb. 18,— DM. 


Tables d’interöts et d’annuites’ Edit&es par le er&dit 
communal de Belgique. 163 S. Bruxelles 1950 (nicht 
im Handel). 


NACHRICHTEN 


Professor Albert Betz 65 Jahre alt 


Am 25. Dezember 1950 beging Prof. Dr. phil. 
Dr.-Ing. e. h. Albert Betz seinen 65. Geburtstag, 

In dem bisher vorliegenden Werk des verdienten 
Gelehrten lassen sich unschwer zwei Hauptrichtungen 
erkennen. Betz nahm ja in der Frühzeit der modernen 
Aerodynamik, im Herbst 1911, als Mitarbeiter von 
Prof. Prandtl seine Tätigkeit in Göttingen auf und 
ließ es sich besonders angelegen sein, die aerodyna- 
mische Versuchstechnik zu entwickeln und zu fördern. 
Die Verbesserung der Windkanäle und ihres Zube- 
hörs, wie der Waagen und der Manometer verdankt 
ihm bedeutende Fortschritte. Nach seinem eigenen 
Zeugnis hat Betz wohl mehr als die Hälfte seiner 
Arbeitszeit auf diese meßtechnischen Fragen ver- 
wenden müssen. Im Mittelpunkt seiner literarischen 
Tätigkeit stehen die Anwendungen der Tragflügel- 
theorie und ihre Fruchtbarmachung für Strömungs- 
maschinen, wie Propeller, Axial- und Radialgebläse. 
In weiten Kreisen ist sein vorzüglicher Artikel über 
„Tragflügel undhydraulischeMaschinen“im VII. Band 
des Handbuches für Physik aus dem Jahre 1927 be- 
kahnt geworden, sowie der spätere Artikel über 
„Applied Airfoil Theory“ im IV. Bande des Durand- 
schen Handbuches ‚‚Aerodynamics Theory‘ aus dem 
Jahre 1935. Aus der letzten Zeit ist von einschlägigen 
Arbeiten die 1948 im Ingenieur-Archiv erschienene 
Abhandlung über den Einfluß der Elastizität der 
Gase auf die Wirkung von Schaufelgittern zu 
nennen; die Strömung durch Flügelgitter pflegt ja 
gern als vereinfachtes Abbild der Vorgänge in Strö- 
mungsmaschinen herangezogen zu werden und daher 
ist die Berücksichtigung des Kompressibilitätsein- 
flusses bei Flügelgittern für die moderne Entwicklung 
besonders dringend geworden. 

Diese wenigen Andeutungen über das reichhaltige 
wissenschaftliche Werk von Betz müssen an dieser 
Stelle genügen. Die Leser dieser Zeitschrift wird sein 
kürzlich bei Springer erschienenes Buch über ‚‚Kon- 
forme Abbildung‘ besonders anziehen, da kein 
anderes zeitgenössisches Werk so wie dieses die prak- 
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W, Weitzel (Prof. a. d. Univ. Kom); Le % Be ch" 
dertheoretischen Physik. 2.Bd.: Struk- 
tur der Materie. XI + 767 S. mit 194 Abb. 
Berlin - Göttingen - Heidelberg 1950. Springer-Verlag. 
Preis geh. 66,— DM, geb. 69,90 DM. 


K..Federhofer (o. Prof. a. d. Techn. Hochschule 
Graz, Dynamik des Bogenträgers und 
Kreisringes. XII+ 179 S. mit 35 Abb. und 
en Harz Wien 1950. Springer-Verlag. Preis ’ 


L. Föppl (o. Prof. a. d. Techn. Hochschule Mün- 
chen) und E, Mönch (a. o. Prof. a. d. Univ. Tucuman), 
Praktische Spannungsoptik VI+ 
162 S. mit 135 Abb. Berlin - Göttingen - Heidelberg 
1950. Springer-Verlag. Preis geb. 21,— DM. 


L, Bieberbach, Einführung in die ana- 
lytische Geometrie. (4. neubearb. Aufl.) 
168 S. mit 43 Abb. Bielefeld 1950. Verlag für Wissen- 
schaft und Fachbuch G.m.b.H. Preis 8,90 DM. 


W,Lietzmann (Prof. a. d. Univ. Göttingen), ‚Wo 
steckt der Fehler?“ Mathematische 
Trugschlüsse und Warnzeichen.183S8. 
mit 120 Fig. Leipzig 1950. B. G. Teubner Verlags- 
gesellschaft. Preis geb. 5,40 DM. 
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tischen Erfordernisse bei der wirklichen Auführung 
konformer Abbildungen betont. 

Seit 1937 war Betz Direktor der Aerodynamischen 
Versuchsanstalt in Göttingen. Nach dem vollstän- 
digen Abbau dieser weltberühmten Anstalt nach dem 
Waffenstillstand fand Betz ein neues fruchtbares 
Wirkungsfeld im Max Planck-Institut für Strömungs- 
forschung in Göttingen. Der Verfasser dieser Zeilen, x 
der an seiner Seite wirken darf, möchte im Namen Bi 
der Fachgenossen dem Jubilar noch viele weitere e 
Jahre freudigen Schaffens wünschen, in denen er 
seine reiche Erfahrung und sein hohes Berufsethos an fe. 
Mitarbeiter und Schüler weitergeben möge. Bu 

Göttingen. - W. Tollmien, 

Freiburg i. Br.: Im Rahmen des wieder auf- 
genommenen Austausches zwischen der Universität 
Freiburg i. Br. und der Technischen Hochschule Karls- 
ruhe hält Prof. Dr.-Ing. Th. Pöschl in diesem x 
Winter eine Gastvorlesung übar „Analytische Mecha- x 
nik“ an der Universität Freiburg. Im Sommersemester N 
1951 wird Prof. Dr.H. Görtler eine Gastvorlesung MW 
en „Grenzschichttheorie“ an I A Eh Karlsruhe 

alten. 


München. Am 6. Oktober 1950 BR im 
76. Lebensjahr der o. Professor für Konstruktions- 
lehre der Dampf- und Verbrennungsmaschinen und. 
Era gR an der Techn. Hochschule NE 

arx 


Es 
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Drückfehlberchi } I: 
W, Lode, ‚Zur thermokinetischen. heo: 
der Elastomere“, Z. angew. Math. Mech, 30 

S. 237/238. Es soll heißen FE 
auf B. 237 Zeile B Yun „außer dem 
ar dem“, 
”’ 11 2 3 Er it J a 

; ta 
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